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We propose a procedure for the determination of the time-dependent ve- 
locity and pressure fields of an unbounded incompressible viscous fluid in an 
external force field induced by an arbitrary number of spheres moving and 
rotating in it as well as the forces and torques exerted by the fluid on the 
particles. Within the completely linearized scheme, we express the velocity 
and pressure fields of the fluid in terms of induced surface force densities and 
derive the explicit form for all quantities contained in these relations not im- 
posing any additional restrictions on the size of particles, distances between 
them, and the frequency range. We show the incorrectness of similar results 
obtained earlier by several authors because these results are expressed in terms 
of nonexistent inverse tensors. We explain the reasons leading to this and pro- 
pose a procedure for the elimination of divergent quantities. In the stationary 
case, using the proposed procedure, we obtained the translational, rotational, 
and coupled friction and mobility tensors for a system containing an arbitrary 
number of spheres up to, respectively, the second, forth, and third orders in 
the dimensionless parameter equal to the ratio of a typical radius of a sphere 
to a typical distance between two spheres. In various particular cases, the re- 
sults obtained in the present paper agree with the well-known results derived 
by other methods. 
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I. INTRODUCTION 



Investigation of hydro dynamic interactions between spheres immersed in an incompress- 
ible viscous fluid is of considerable interest for a wide class of problems of physics of suspen- 
sions and colloidal crystals. Due to the long-range nature of the hydrodynamic interactions, 
their account is essentially important for the study of these systems, for which many-particle 
interactions should be taken into account. As a rule, for the study of these interactions, the 
linearized Navier-Stokes equation for the fluid is used ||T|. Usually, the problem is reduced 
to the determination of the forces and torques exerted by the fluid on particles moving and 
rotating in it with given translational and rotational velocities. According to the classical 
approach used for one particle 0-^, to find these forces and torques, first, it is necessary 
to determine the velocity and pressure fields of the fluid induced by particles moving and 
rotating in it. However, in the case of several particles, this problem is extremely compli- 
cated. The corresponding exact relations for the forces and torques exerted on particles due 
to hydrodynamic interactions between them were derived only for the particular stationary 
case of two spheres moving along the line passing through their centers with equal or dif- 
ferent velocities |^ and perpendicular to this line and rotating along the axis perpendicular 
to the direction of motion of the spheres and the line connecting their centers P,p!0[ . 

For this reason, several methods aimed at the determination of approximate solutions 
have been developed. Among these methods, there is the well-known classical method of 
reflections first proposed by Smoluchowski for analysis of the forces exerted by the fluid on n 
spheres moving in it with constant velocities. Later, this method was used for the solution of 
many stationary problems of hydrodynamic interactions of particles (mainly, two particles) 
in an incompressible viscous fluid (the detailed review of results is given in [jl|) including 
also the cases of permeable spheres [p!T| -[T3[| and mixed slip-stick boundary conditions at the 
surfaces of the spheres [p!4|-[T8|| . The corresponding results are presented in the form of power 
series in the dimensionless parameter a equal to the ratio of a typical radius of a sphere to 
a typical distance between two spheres calculated to a certain order of this parameter. 
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Since the method of reflections and its further modiflcations are based on the solution 
of the corresponding boundary-value problems, their direct application to the solution of 
nonstationary many-particle problems seems to be rather problematical. Moreover, even 
in the stationary case, the procedure of determination of the velocity fleld of the fluid 
induced by particles becomes essentially more complicated with increase in the number 
of particles. For this reason, the methods of reflections are usually used for two spheres 
in the stationary case. In Mazur and Bedeaux proposed a new method (called the 



method of induced forces) for the determination of the force exerted by the fluid on a 
single sphere moving with time-dependent velocity in the case where the fluid moves with 
nonstationary and nonhomogeneous velocity. Later, this method was developed for the 
determination of the friction pO[ and mobility tensors in the case of the stationary fluid 



containing an arbitrary number of particles, furthermore, the results for the mobility tensors 
were generalized to the nonstationary case p2[. At flnite frequencies, the expressions for 



the mobility tensors are obtained up to the third order in two dimensionless parameters, 
namely, the parameter a, which is typical of the stationary case, and the parameter k, 
which is proportional to the ratio of a typical radius of a sphere to the penetration depth 
of transverse waves. The method of induced forces essentially differs from the method of 
reflections because it is not based on the necessity of the knowledge of the explicit form for 
the fluid velocity induced by particles. 

Other methods for the solution of problems of hydrodynamic interactions between par- 
ticles in the fluid were proposed in [p^ -p9| . Despite the fact that the solution of nonstation- 
ary many-particle problems in [^,0 is also based on the introduction of unknown induced 
forces, this approach essentially differs from the method of induced forces used in p0|-p2 



Indeed, within the framework of this approach, the forces and torques exerted by the fluid 
on particles as well as the velocity fleld of the fluid induced by these particles are expressed 
in terms of the introduced induced surface forces, the explicit form for which was determined 
up to the third order of two dimensionless parameters mentioned above, while the method 
of induced forces in p0|-[2^ is developed in such a way that the mobility and friction tensors 



(not the velocity field of the fluid that cannot be found with the use of this method) are 
found without determination of the explicit form of the introduced induced surface forces. 
However, the final results presented in p5| , |26[] are expressed in terms of multiindex tensors, 
the explicit expressions for which are not given and that are only defined as tensors inverse 
to certain rather complicated multiindex tensors. For this reason, the question of the agree- 



ment between the results obtained in ||25| , |26[| and is open. Furthermore, this problem 



remains unsolved even for the stationary case. Even in the simplest case of a single sphere 
rotating with constant angular velocity in a fiuid, the fiuid velocity followed from the general 



relations given in p5| , ^ is also expressed in terms of the inverse tensor and, hence, cannot 
be reduced to the classical result for the fiuid velocity induced by a rotating sphere [jl|,|^ 
without the knowledge of the explicit form for the inverse tensor. 

The aim of the present paper is to develop a method for the solution of various problems 
(both stationary and nonstationary) of hydrodynamic interactions between any number of 
particles immersed in an unbounded incompressible viscous fiuid in an arbitrary force field 
such that, as opposed to [^,^, the results obtained by this method (the velocity and 
pressure fields of the fiuid induced by these particles, the forces and torques exerted by the 
fiuid on the particles, various mobility tensors, etc.) coincide with the corresponding results 
obtained in particular cases by other methods as well as to explain the reasons that do not 
enable one to represent the results given in in the conventional form. 



In Sec. 2, we reduce the problem of n spheres in an unbounded incompressible viscous 
fiuid in an arbitrary force field to an equivalent problem for this fiuid without particles in 
an efficient force field. Unlike the generally excepted procedure of the zeroth extension of 



the external force |[T9|j20|j23| --p6| and the fiuid pressure 0-|2^ to the domains occupied by 
the particles, we use another extensions for these quantities and discuss possible problems 
connected with different methods of extensions of the required quantities. We give relations 
for the total forces and torques exerted by the fluid and force flelds on particles. 

In Sec. 3, within the framework of the completely linearized scheme (both with respect 
to the fluid velocity and the velocities of the particles), we express the required distributions 



of the fluid velocity and pressure in terms of the unknown induced surface forces distributed 
over the surfaces of the particles. The results are obtained without imposing any additional 
restrictions on the size of particles, distances between them, and frequency range. All 
required quantities contained in the relations for the velocity and pressure flelds of the fluid 
are given in the explicit form in terms of special functions of dimensionless parameters. We 
derive the system of algebraic linear equations in the unknown harmonics of the induced 
surface force densities in Sec. 4. 

In Sec. 5, we consider the stationary case. Taking into account the obtained explicit 
form of the quantities in the system of algebraic equations, which is similar to the system 



given in ||3^, we solve this system by the method of successive approximations using the 
system of noninteracting particles as a zero approximation. In this approximation, we show 
that the determinant of the system of equations corresponding to the harmonics with / = 1 
is equal to zero. Therefore, the inverse tensor, in terms of which the relations for the forces 



and torques exerted by the fluid on the particles given in |3y] for the stationary case as 
well as analogous results and the velocity fleld of the fluid induced by the particles given in 
for the nonstationary case, does not exist. We analyze the reasons for this fact and 



show that this means that the required induced surface forces can be determined only up to 
arbitrary potential components that have no influence on the fluid velocity. We formulate 
the additional conditions that enable one to uniquely determine the induced surface force 
densities. Using the proposed procedure, we obtained the translational, rotational, and cou- 
pled friction and mobility tensors for a system containing an arbitrary number of spheres up 
to the second, forth, and third orders in the dimensionless parameter o", respectively, which 
in various particular cases, agree with the well-known results obtained by other methods. 
Within the framework of the considered approach, we formulate an algorithm for the de- 
termination of the velocity and pressure flelds of the fluid as well as the forces and torques 
exerted by the fluid on particles corresponding to a given power of the parameter a. The 
details of calculation of several integrals of the products of three and two Bessel functions 
necessary for determination of the hydrodynamic interaction tensors are given in Appendix. 



II. GENERAL RELATIONS 



We consider an unbounded incompressible viscous fluid in a certain external force field. 
In the linear approximation, the fluid is described by the linearized Navier-Stokes equation 

P^^^ + div P{r, t) = F-*(r, t) (2.1) 

and the continuity equation 

dwv{r,t)=0. (2.2) 

Here, P(r, t) is the stress tensor of the fluid with the components 

P,j{r,t) = 5,,p{r,t) - r, (^^^ + ^^^) > = x,y,z, (2.3) 

p{r, t) and v{r, t) are the hydrostatic pressure and velocity fields of the fluid, p and r] are its 
density and viscosity, respectively, Sij is the Kronecker symbol, and F^^^{r, t) is the external 
force acting on a unit volume of the fluid. 

We represent the external force F'^^*{r, t) as a superposition of the potential F^^'^^^{r, t) 
and solenoidal F^*°'^^'^*(r, i) components 

F^^*(r, t) = F(f)"^*(r, t) + F(^°')^^*(r, t), (2.4) 
F^>'''{r,t)^-pVip{r,t), (2.5) 

where ip{r, t) is the potential of the external conservative force. 

Assuming that the quantities in Eqs. (2.1)-(2.3) can be expanded in the Fourier integral 
of the form 
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A{r, t) = — ^ Jduj Jdk A{k, cu) exp(i(fc • r - cut)), (2.6) 



-oo — oo 



(27r)^ 

we represent the frequency Fourier transform of the solution of Eqs. (2.1)-(2.3) as follows: 

(Zn]-^r} J K'^ -\- K,-^ 



— OO 

p(r, ^) = -7^ Jdk !^P(^ k . F-\k, u;) + p-f{u;) 



{2nf _ 

-pip{r,u;) + p'''f{u;), (2.8) 



where 

^isol)ext^f^^ CJ) = (/ - UkUk) ■ F^^\k, u) (2.9) 

/ is the unit tensor, rik = k/k is the unit vector directed along the vector k, hi = ^Juj /{iu) = 



(1 — isignuj)/6, 6 = y2z//|c<j| is the depth of penetration of a plane transverse wave of 
frequency uj created by an oscillating solid surface into the fluid u = t]/ p is the kinematic 
viscosity of the fluid, v'^^^^cu) = 2TT5{ij)v^^^ and p^"^'^ {uj) = 27r5(u;)p*"-^, where 5{iu) is the 
Dirac delta-function, i;*"-^ and p*"-^ are certain constant velocity and pressure at infinity 
defined by the condition for the fluid at infinity where the external force field is absent. 
According to relation (2.7), the fluid velocity is determined only by the solenoidal component 
of the external force, which is the natural consequence of the continuity equation (2.2). The 
fluid pressure is determined by the potential component of the external force. 

Let, in the considered fluid, homogeneous macroscopic spheres of radii Oq and masses 
rria, where a = 1, 2, . . . , iV, be present. At the time t, the position of the center of sphere a 
is defined by the radius vector Ra(t) relative to the fixed Cartesian coordinate system with 
origin at the point O (in what follows, the system O). In parallel with this fixed system, we 
also introduce A^ local moving Cartesian coordinate systems with origins Oa at the centers 
of spheres coinciding with their centers of mass (in what follows, systems Oa) so that sphere 
a does not move relative to system O^- Therefore, the radius vector r of any point of the 
space can be represented in the form r = Ra + r^, where is the radius vector of this point 
relative to the local system Oa- The motion of the spheres is described by the equations 

m.^^ = Fi^\t), (2.10) 

/a^ = T-(t), (2.11) 

where 

_ dRajt) 

dt 

and 0,a(t) are, respectively, the translational and angular velocities of sphere a, la = 
{2/5)maaa is its moment of inertia, F^°^{t) is the total force and T^a^(t) is the total torque 



(here and in what follows, all torques corresponding to particle a are considered relative to 
its center) acting on sphere a 

F'f{t) = Fl^\t)+Fi{t), (2.12) 

r:\t) = Ti^\t) + Ti{t), (2.13) 

where F'^\t) and T'^^(t) are the force and torque acting on sphere a due to the external 
force field and F{^{t) and T{(t) are the force and torque exerted by the fluid on sphere a 



Fiit) = - J P{r,t)-n^dS^, (2.14) 
Tiit) = - J {{r-R^{t))xP{r,t))-n^dSo., (2.15) 

Sa{t) is the surface of sphere a at the time t and na is the outward unit vector to this 
surface. 

The fluid is described by the same equations (2.1)-(2.3) as in the case of the absence 
of particles but defined only in the domain outside the spheres Tq, > Oa, a = 1, 2 . . . , iV, 
where = \ra\ = |r — Ra(t)\. 

The problem is to determine the velocity and pressure fields of the fluid induced by 
moving particles in it as well as the forces and torques exerted by the fluid on the particles. 
For the solution of this problem, we use the method of induced forces proposed in |[T^ and 



developed in p0|-[2^. We extend the quantities p{r,t), v{r,t), and F^^^{r,t) defined in 
Eqs. (2.1) and (2.2) for > aa to the domains Tq, < Oq, occupied by the spheres as follows: 

vir,t) = Uair,t) = Uait) + in^it)xr^), < a^, (2.16) 
p(r,t) = -py,(r,t)+p-^-p(^r.^^j, r„ < a„. (2.17) 

We extend the external force field F^^^{r,t) given for > to the domains < 
so that in the entire space, it is described by the same analytic expression as in the domain 

> Oq. (For short, we call this extension an analytic extension.) The nonzero extension 
of the fluid pressure defined by (2.17) and the analytic extension of the external force field 
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differ from the usually accepted zeroth extension of both the fluid pressure [pJl-p^ and the 
external force field p3|-p^ to the domains < a^. It is worth noting that the considered 
problem can be solved for various extensions of the quantities p{r,t) and F^^*{r,t) to the 



domains < including the zeroth extension |p!9|-[26|] and the extensions used in the 
present paper. However, the complexity of calculations and the possibility of interpretation 
of certain results essentially depend on the specific method of extension. 

If the fluid pressure is extended to the domains < according to relation (2.17), then 
the function p{r,t) defined in the entire space has a discontinuity at the surfaces = Oq. 
For this reason, in what follows, we write the fluid pressure at the point r = + cIq as 
p{Ra + aa + 0,t). 

Note that for the analytic extension of the external force field F^^^{r,t) to the domains 
Tq, < Qa, the first two terms in (2.17) represent the pressure of the unbounded fluid without 
particles in the external force field F^^^(r,t). In the case where the fluid velocity satisfies 
the stick boundary conditions at the surfaces of the spheres |]1]^J§] (this case is considered 
in the present paper), condition (2.16) ensures the continuity of the function v{r,t) given 
in the entire space at = a^. 



In P^-pB[, the external force field F^^ (r, t) is extended to the domains < by zero 



and the stress tensor P{r, t) in these domains is defined as follows: 

divP(r,t) = -p^^^^^, r«<a„. (2.18) 

We note that in the case where the fluid velocity v{r,t) is extended to the domains 
fa < O'a according to relation (2.16) (in fact, this is used in p^-p6[), the stress tensor 
P{r,t) for Ta < da must have another analytic representation than (2.3). Otherwise, if the 
stress tensor P{r,t) in the domains Tq, < Oq, is defined by relation (2.3), where the quantity 
v{r,t) is defined by relation (2.16) in these domains, then for < a^, we have 

P{r,t)=p{r,t)I, (2.19) 

where the quantity p{r, t) for r„ < depends on the method of extension of the fluid 
pressure p(r, t) defined for > to < a^. However, for any extension of the fluid 



pressure p{r,t) to the domains Tq, < a^, according to (2.19), the left-hand side of Eq. (2.18) 
contains the potential vector Vp(r,t), while the right-hand side of this equations contains 
the solenoidal vector -p^ ((f2a(i) x r)) because 

{n^{t) X ^) = ^ rot((rj^(t) X r) X r). (2.20) 

This means that the extension of the fluid velocity v{r, t) to the domains < according 
to (2.16) eliminates the possibility of extension (2.18) (used in [p^ -|26[|) for the stress tensor 



P{r,t) defined by relations (2.3) in the entire space to these domains for the case of the 
time-dependent angular velocity flait). 

The quantities v{r, t) and p{r, t) extended to Tq, < Qq, according to (2.16) and (2.17) and 
F^^^{r,t) analytically extended to these domains are given in the entire space and satisfy 
the continuity equation (2.2) and the equation 

^dv{^ + div P(r, t) = F^^\r, t) + F™'^(r, t) (2.21) 

given in the entire space. Equation (2.21) differs from Eq. (2.1) by the additional term (the 
induced force density) F^^'^{r,t) that appears due to the extension of the linearized Navier- 
Stokes equation given for > to the domains < aa- For the analytic extension 
of the external force field, the quantity F'^^^{r,t) in Eqs. (2.1) and (2.21) is described by 
the same analytic expressions, while, for the zeroth extension of the external force field to 
Va < tta ||T9| , [23| -p6| , the analytic expressions describing the quantity F^^^{r, t) in the domains 
Va > aa and Tq, < are different. Therefore, in the latter case, the procedure where the 
quantity i^*"°'(r,t) is formally considered to be equal to zero in Eq. (2.21) is not equivalent 
to the removal of the particles from the fluid because Eq. (2.1) given in the entire space and 
Eq. (2.21) with F^"^(r, t) = describe the unbounded fluid in various force fields, namely, in 

N 

F^^ (r,t) and in J2 'S>{ra — aa)F^^ (^;^); where 6(x) is the Heaviside function, respectively. 
The induced force density F^'^'^(r,t) is presented in the form 

N 

F'^^^ir, t) = J2 K\r, t), (2.22) 
«=i 
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where F^^'^{r,t) is the induced force density for particle a. This quantity has the volume 
j7i(y)mdj^^^ t) and surface F^^'^"''^{r, t) components given, respectively, inside the volume \4 = 
(4/3)7raj^ occupied by the particle a and on its surface Saif) 



F'^\r, t) = Fp"'^(r, t) + Fp"'^(r, t). (2.23) 

Here, 

= e(a, - {n^{t) x - F(^°')^^*(r, t)| , (2.24) 

Fp"'^(r,t) = Jdn^dir - R^it) - a„)/Ja„,t), (2.25) 

where = (oq,, 9a, ^Pa) is the vector directed from the center of sphere a to a point on its 
surface characterized by the polar 6a and azimuth angles in the local spherical coordinate 
system Oa, dQa = smOadOadipa is the solid angle, and fa{o,a,t) is a certain unknown 
density of the induced surface force distributed over the surface Sa(t). 

Unlike ||21|, the volume component i^^^^*"'^(r, t) of the induced force does not contain 
the translational velocity Ua{t) of particle a, which is caused by the different extensions 
of the fluid pressure to Tq, < aa used in the present paper and in In addition, owing 
to extension (2.17) for the fluid pressure to Tq, < a^, i^^^-'*"''^(r, t) is defined only by the 
solenoidal component of the analytically extended external force field. Thus, in the particular 
case of conservative external forces, i^^^^*"°'(r, t) is independent of external forces. 

It is easy to verify that the representation of the induced force density in the form 
(2.22)-(2.25) as well as relations (2.16) and (2.17) and the analytic extension of F''^*(r,t) 
to Ta < tta insure the validity of Eq. (2.21) in the entire space. In this case, 

div P(r, t) = -P^^^ + F^P^'^\r, t) + Fp"'^(r, t), < a„. (2.26) 

Using relation (2.26), we can represent force (2.14) and torque (2.15) exerted by the fluid 
on sphere a as follows: 

Fiit) = Fait) + F^:{t) - F^^^-'^'it), (2.27) 
Ti{t) = Ta{t)~f^:^'''\t), (2.28) 
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where Fait) and Taif) are, respectively, the force and the torque exerted on the sphere a 
due to the induced force density distributed over its surface 

F^it) = - jdrF'-i^'^^r.t) = - JdQ^ f,{a^,t), (2.29) 
T^{t) = - jdr (r, X FP"^(r,t)) = - jdQ^ (a, x /Ja,,t)), (2.30) 

F^^ (t) and T^f ^ (t) are, respectively, the force and the torque acting on the fluid sphere 
occupying the volume Va instead of spherical particle a due to the potential component 
F^^^^^^{r,t) of the external force field analytically extended to this domain 

= jdrF^P^^^\r,t), (2.31) 

f^^''\t) = Jdr (r„ X F(^)-*(r, t)) , (2.32) 

and 

FTit) = (2-33) 

is the inertial force that is necessary to be applied to a fluid sphere of volume Va in order 
that it move with the acceleration dUa{t)/dt |T9|, where rha = pV^ is the mass of the fluid 
displaced by particle a. 

In the particular case of the time-independent homogeneous gravitational force field 

ip{r, t) = ip{r) = -gr, F{^oi)ext^^^ _ ^2.34) 

we have 



(2.35) 

TTit) = 0, T:'-{t) = 0, 



(t) = Ff = m^g, (t) = -Fi 



where F^ = —rhag is the Archimedes force acting on particle a and g is the acceleration 
of gravity. Therefore, in the considered case, F'^^{t) — F^ (t) = {nia — rha) g in relation 
(2.12) is the gravity force of particle a corrected for the buoyancy force. 

In view of relations (2.29) and (2.30), the forces and torques exerted by the fluid on 
particles are completely determined by the surface induced force densities. If it is necessary 
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to determine only the forces Fa{t) and the torques Ta{t) in the hnear approximation with 
respect to the velocities of the fluid and the spheres, it is convenient to use the procedure 
proposed in pO|- p^] , which does not require the explicit form of F^^^^'^'^{r,t). On the basis 
of this procedure, for an arbitrary number of spheres in a fluid, one obtains the mobility 
tensors of particles that move and rotate in the fluid both in the stationary case pD|pT| 
and with regard for the time dependence ||2^. In a more general case where the velocity 



field of the fluid should be determined, the problem becomes essentially more complicated 
because the explicit form of the induced surface force densities must be obtained. In the 
nonstationary case, this problem is considered in ||25| , |26|1 where the zeroth extension of the 



external forces to the domains occupied by spheres is used and the stress tensor in these 
domains is defined by relation (2.18). 

III. REPRESENTATION OF THE REQUIRED QUANTITIES IN TERMS OF 
HARMONICS OF INDUCED SURFACE FORCE DENSITIES 

To determine the velocity and pressure fields of the fluid in the presence of spheres in 
it, we note the same structure of Eqs. (2.1) and (2.21). Therefore, we can use solution (2.7), 
(2.8) of Eqs. (2.1)-(2.3) with the substitution F"^'*(fc,cj) + F'''\k,uj) for F"^'*(fc, cu). 

We consider the problem linearized both with respect to the velocity of the fluid and the 
velocities of the spheres. Within the framework of this approximation, we neglect the time 
dependence of the positions of centers of the spheres defined by Rait) and their surfaces 
Sa{t). This means that the spheres do not displace considerably for considered time intervals 
(for details, see |22|). This enables us to represent F^f )*"°'(fc, cu) and F(,^)™^(fc, cu) as follows: 



FP"'^(fc,cu) = Jdnaexp{-tk-{Ra + aa))fMa,uj), (3.1) 
F^^^-'{k,co) = -!^l^X^^^{^a{co) X fc) +Fi^°'^^"(/c,o;)|exp(-zfc.i?.), (3.2) 

where C,a = Qirrjaa is the Stokes friction coefficient for a sphere of radius a uniformly moving 
in the fluid, ba = nda is the dimensionless parameter that characterizes the ratio of the 
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radius a of sphere a to the length 5 [\ba\ — V2aa/S)], jn{x) is the spherical Bessel function 
of order n, and 

F^r'^""\k,u;) = Jdr^ exp {-tk ■ r«) Q{a^ - r«) F^'"'>^\R^ + r„,u;). (3.3) 

Note that the quantity ik,uj) defined by relation (3.3) differs from the Fourier 

transform F^^''^'^^^(k,uj) of the solenoidal component of the external force field because the 
integral in (3.3) is taken over the finite space domain Tq < Cq,, while for F^*°'-*'^^*(fe, a;), the 
Fourier integral is taken over the entire space. 

By using relations (2.7), (2.8), and (3.1)-(3.3), we can represent the required Fourier 
transforms for the quantities v{r,u!) and p(r,a;) defined in the entire space as follows: 

v{r,u)^v^^\r,u) + v''''^{r,u), (3.4) 
p{r, uj) = (r, u) + p'^^^ir, cu). (3.5) 

Here, the quantities p^'^\r,u!) and v^'^\r,cj), where 

t;(°)(r,a;) = t;(°)^°'(r,cu) + 'y*"^(cu), (3.6) 

— oo 

defined at any point of the space, for Tq, > a^, are, respectively, the pressure and the velocity 
of the unbounded fluid in the absence of particles deflned by relations (2.8) and (3.6)-(3.7). 
The quantities i;*"'^(r,a;) and p^'^'^ir^uj) can be represented in the form 

N 

v'^\r,uj) = Y.v'^\r,uj), (3.8) 

13=1 

N 

p^^\r,u) = Y.p^^\r,u), (3.9) 

/3=1 



where 



<(r-,u;) = v^P'-\r,u) + vf'''\r,u^l (3.10) 
p^^\r,i.)=pf^'''\r,i.) +pf'''\r,i.). (3.11) 
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For Tq > Oq,, v'"''^(r, o;) and p*"'^(r,a;) are, respectively, the fluid velocity and pressure 
induced by all particles. Here, 



— oo 

(3.12) 

pf^''"^{r,uj) = - Jdk ^^P(^^ — Jdfi^ exp(-ifca^) A; • fp{af),uj). (3.13) 

— oo 

Relations (3.12) and (3.13) are valid for any r. For > a^, the quantities v^^^*"^'^(r, a;) and 
^{S)md^^^ o;) can be interpreted, respectively, as the fluid velocity and pressure at the point 
r generated by the induced surface force F^^^'^''' {r , lj) distributed over the surface of sphere 



The quantity 

S{k,u)^-^^^^{I-n,n,) (3.14) 

is the dynamic Oseen tensor. 

The quantities 'y^^^*"'^(r, a;) and p^^^^^^ {r , ou) defined in the entire space, in the domain 
Tq, > tta, are, respectively, the fluid velocity and pressure at the point r generated by the 
induced volume force l^^^^*"°'(r", a;) given in the volume V/j occupied by particle (3, moreover, 

* J*!3i(!^Jpl^fc.F<-)-(fe,^). (3.15) 

— oo 

Thus, the contribution to the fluid pressure due to the volume induced forces is nonzero 
only if the considered system is in a nonconcervative force field. 



It is convenient to represent the quantity t;L^'**"'^(r, a;) as the sum of two terms 



v^P'''\r,uj) = vf'''\r,uj) + ^;^^°'^*""(r,^), (3.16) 



where the first term is caused by the rotation of particle /5 



(T)ind, X 22 1 7 exp(ik-(r-R/3))j2(Kap)f ^ 
vy (r,^) = -^^pb,--^^ jdk ^(^Pi^) X k) (3.17) 



-oo 



'3ri {2nY 

and the second is caused by the solenoidal component of the external force filed 
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n ^'''^^ = "{2^ J WT^^ (/-nfcnfc)-F^ {k,uj). (3.18) 



m 



Note that for the zeroth extension of the external forces to the domains Tq, < aa used 
20| , ^3| -p^ ] , the interpretation of v^'^\r,iu) as the nonperturbed velocity of the fluid, i.e., 



the fluid velocity in the absence of the spheres, is not a quite correct because, in this case, 
v^^\r,uj) is the velocity of the fluid not containing particles but in the external force field 
J2 ©('"a ~ O'ajF'^^ {f^'t) instead of F^^ {i^,t) given in the entire space. The same remarks 
are also valid for the interpretation of p^^\r,Lj). 

We expand the induced surface force densities fpi^ap.uj) in the series in the spherical 
harmonics 

//3(«/3, = lX= Uj)Ylm{0p, ^p), (3.19) 

^V^ Im 

oo I 

where, for short, the symbol means X) S and the arguments a^g in f f^^ij^{ap,uj) is 

Im 1=0 m=—l 

omitted in what follows, and the spherical harmonics Yjm(^)V^) are defined as follows |pi|| : 
Yi„,{e, if) = ^^^Pi^0,^(cos^), / = 0, 1, . . . , -/ < m < /, 



jTi — I m I 

QiM = (-1) 



!i±i«pH(.). (3,20) 

\ 2 [l+\m\)\ 

and P[^{x) is the associated Legendre polynomial. 

By using expansion (3.19) and relations (2.29) and (2.30), we can represent the Fourier 
transforms of the force Fa{uj) and the torque Ta{uj) exerted by the fluid on particle a in 
terms of Fourier harmonics of the induced surface force density fdmi^) as follows: 



rn=-l 

where 



Faioo) = -/.,oo(^), (3.21) 
Ta{u;) = E {^m X /„,i^(^)) , (3.22) 



Co = e,, e±i = [isy ± e^.) , (3.23) 
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i.e., up to the factor (—1)"*, where m = 0, ±1, are the cychc covariant unit vectors [0, and 



62;, and are the Cartesian unit vectors. 

Thus, the problem of determination of the forces and torques exerted by the fluid on 
particles is reduced to the determination of only the zeroth and first (with respect to /) 
harmonics of the induced surface force densities [ P^ , P^ ,Pm . 

For any point of the space, its position defined by the radius vector r relative to the 
fixed system O can be represented in the form r = + Va, where Ra is the radius 
vector that defines the position of the center of certain sphere a and Va is the radius vector 
directed from this center to the point at hand. Since there are particles, there exist 
different representations for r. Among these representations, in the case where the point 
of observation belongs to the domain occupied by the fluid, it is convenient to take r„ 
corresponding to the minimum difference rp — ap, where (3 = 1,2, . . . , N, which means 
the minimum distance from the point at hand to the surface of sphere a. If the point of 
observation lies inside the domain occupied by sphere a, then is the radius vector directed 
from the center of this sphere to this point. 

In this case, for Tq, > a^, the quantities v^p^'^'^{Ra + ra,uj), p^P^^'^{Ra + ra,uj) and 
Vp^^^"''^{Ra + Va, 00), p^^'**"°'(-Rq, + Va, uj) , where f3 = 1,2, . . . , N, are the velocities and the 
pressures of the fluid at the point r = Ra + induced, respectively, by the surface and 
volume forces distributed over the surface and inside the volume of particle (3 with the center 
defined by the vector i?^. The particular case P = a corresponds to the velocities and the 
pressures of the fluid generated at this point, respectively, by the surface and volume forces 
distributed over the surface and inside the volume of the particle that is the most closely 
situated to this point. This enables us to separate contributions to the induced velocity and 
pressure of the fluid caused by the closest particle and more distant ones. 

At any point r of the space, we can represent the quantities defined by relations (3.12), 
(3.13), (3.15), (3.17), and (3.18) as expansions in the spherical harmonics 

^^S>nd^^^ + u) = 2y/7C ^ vf]^'^{Ra, r„, Uj)Yim{Oa, fa), (3.24) 

Im 
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p(S)^nd^j^^ + Va, Uj) = 2^/^ Y.PfJm'^i^cy^ r^, Uj)Yirni9a, fa), (3-25) 

Im 

v''^^''''^{Ra + ra,uj) = 2y/n^v^^}l^'^{Ra,ra,u;)Yim{0a,(Pa), (3.26) 

Im 

^(so«>nd^^^ + ra,u}) = 2x/7r ^ v5j7m"'^(-^«,ra,a;)Y'i,„(6'a,93a), (3.27) 

Im 

pip^nd^j^^ + Va, uj) = 2V7r ^ p^^)^'"^(Ha, r^, a;)Yi^(ea, V'a), (3.28) 

Im 

where Tq, = (ra, 9a, (Pa) is the radius vector from the center of the particle closets to the 
considered point if > or the radius vector from the center of the particle to the point 
of observation lying inside this particle if Tq, < Ga and 9^ and are, respectively, the polar 
and azimuth angles of this vector in the local spherical coordinate system O^- The expansion 
coefficients in series (3.24)-(3.28) of the corresponding quantities in the spherical harmonics 
Yim{9a, fa) Can be presented in the following form: For (3 ^ a, we have 

vf,lZ'iRo.,ra,uj) = E T^ai^a,ap,Rap,uj) ■ fp,i,mM^ (3-29) 
pfJ:i'{Ra,ra,u;) = E D\lZ{ra,afs,Rap) ■ /;3,/.^.(^), (3-30) 

l2m2 

4!tm.iRa,ra,u;) = l^pblJ:{-iy^Pi,2,i2ira,ap,Rap,u;) (fi^^) x ^^,00) 

l2m2 

XYi:^^{Qap,^af3), (3.31) 

u;) = -2V^ jdrpQ{ap - r^) ^ T^^ (r., r^, R^p, u) 

■ F(-')-*(i2^ + r^,^)l^:^^(^^,(^^), (3.32) 
pS2^(J^a, ^a, ^) = 2^/^ /c^r-/3e(a^ - Tp) Y: D\\Z\{ra, rp, Rap) 

12^2 

. F(-')-^(i?^ + r^,u;)r,:„^(e^,(^^), (3.33) 

where Ra/s — Ra — R(3 = (-Ra/3, ^ap, ^ap) IS the vector between the centers of spheres a and 
(5 pointing from sphere (5 to sphere a in the spherical coordinate system O, 

T\lZ\{ra,rp,Rap,^) = Y.T\Z\Mi^o.,rp,Rap,u^)YUQap,^ap), (3.34) 

= Fi,i,,i{r^,rp,Rap,uj)K\lZlim, (3.35) 
27 A;2 

Fi^i^^i(ra,rp,Rap,u;) ^ — Idk ——^ ji^{kro)ji^{krp)ji{kRap), (3.36) 
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T^l2m2 _ Ji-h+l 
limijlm 



jdVtk {I - rikUk) Yi*^^ {9k, ipk)Yi^m2iOk, ifk)yL{Ok, ^k), (3.37) 



Im 



limi ,lm ' 



12^2 
l\m\,lm'> 



IT 



Ci^i^,i{ra,rp,Rap) = - dk k ji^{kra)ji2{krp)ji{kRap), 



(3.38) 
(3.39) 
(3.40) 



rl2m2 



^^limi,lm ^ 





h-h+l-l 

2 7,. k 



jdVtk rik Yilra^iOk, (pk)Yi^m2{0k, (pk)Yi*^{ek, (fik), (3.41) 



Z f K 

Phi,, lira, a Rap, uj) = — dk -——^ ji^{kra)ji2{kap)ji{kRap). 

' 



(3.42) 



hmuim ^^'^ hmulm defined by relations (3.37) and (3.41) can be 



The quantities K\^^'^ ,^ and w'''^"^^ 
represented in the exphcit form in terms of the Wigner 3j-symbols |22|,^. (For the tensor 
^hml imi ^^^^ representation is given in ||3^.) The corresponding relations have the form 



T^l2m2 _ 



where 



V'(2/i + l)(2/2 + l)(2/ + l) (-1)™!+™ 

/ 

V 







V 



-nil m 



k= 2 j — jmin 



2 / L 








X 



2 / L 

k —m m — k 



—nil m2 k — m 



J ) 



k,k,l> 0, 



L = l + 2j, 



I, if / = 
0, if 1 = 1 
4, if l>2, 



(3.43) 



Ko = ^- (-e^e^ - ByBy + 2e^e^) = y 3 (eie_i + e_iei + 2eoeo) 
Ki = e^e^ + e^e^ - i [sye^ + e^Cy) = -a/2 (e_ieo + eoe_i) , 
K I = -KI = -72 (eico + eoei) , 

K_2 = -ft:; = 2eiei, 



(3.44) 
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and 



(imi ,lm 



X 



.^(2/i + l)(2/2 + l)(2/ + l) (-1) 



( 1 I l\ ( h h l \ ( 1 











/ L 
k —m m — k 



k= 1 J — Jmin 



V 



-nil ^2 k — m 



(3.45) 



where 



l + 2j- 1, 



1, if / = 
0, if / > 1. 



Relation (3.43) coincides with relation (28) given in up to the factor 

(_l)(™i-|m-il+"^2-|m2|+m.-|m|)/2^ which is caused by the different definition of the spherical 



harmonics Yi.in(9, (p) for m < in the present paper and in [30 
In view of the properties of 3j-symbols, we have 



K\Zlim^(^ only for I = I, + I, - 2p, / > 0, 



(3.46) 



where p = -1,0, 1, . . .,Pmax, Pmax = niin([(/i + l2)/2] , 1 + min(/i, Zg)), [a] is the integer part 
of a, and min(a, b) means the smallest quantity of a and b, and 



W\IZIm^^ only for / = + _ 2p + 1, / > 0, 



(3.47) 

where p = 0, 1, . . . ,Pmax, Pmax = mm[[{li + I2 + l)/2] , 1 + min(Zi,/2))- 

According to (3.46) and (3.47), the sums over / in relations (3.34) and (3.38) contain 
only terms with I satisfying these conditions. Therefore, it is necessary to investigate the 
quantities Fi-^i^^i{ra, rp, Rap, uj) and Ci^i^^i{ra, rp, Rap) only for the values of / given by (3.46) 
and (3.47), respectively. According to relations (3.29), (3.30) and (3.32), (3.33), the corre- 
sponding quantities in them are defined by general relations (3.34)-(3.41) for rp = ap and 
< ^'/s < a^, respectively. 

In the particular case / = 0, the general relations (3.43) and (3.45) for the quantities 
K\lZ,im and Wll'^li^ are simplified to the form 
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T^l2rn2 _ 



dVtk {I - UkUk) Yil^^ {0k, (pk)Yi2m2 ^k) 



3v^ 







k=-2 



—mi 1712 k 



(3.48) 



and 



^fe.oo = JdnknkY*^^{ek,^k)Yi,m2{0k,'Pk) 



j^h—l2~l 



20F 



(2Zi + l)(2Z2 + 2)(-l) 



mi 



^1 ^2 1 ^ 







k=-l 



(3.49) 



^ —mi 777.2 ^ J 

Taking the properties of 3j-symbols into account, we obtain that 



J^l2m2 



Zimi,00 



^ only for ^2 = + 2p > 0, where p 



0, 1, if Zi = 0, 1 
0,±1, if/i>2, 



(3.50) 



and 



"^hZfio 7^ only for ^2 = + 2p + 1 > 0, where p = < 



0, if /i = 0, 1 



(3.51) 



0,-1, ifZi>l. 

According to (3.51), the sum over I2 in relation (3.31) contains only terms for I2 = 
Zi ± 1 > 0. This means that the quantity Pii2,i2{fa-iCip,Rap-iUj) contained in this relation 
must be determined only for these values of 12- 

In the domains < Rajj — Tjj and > Rap + ^'p (far from both particles), where, 
according to (3.32) and (3.33), rp < ap, the integrals in (3.36) and (3.40) that define 
the quantities Fi^i^^i{ra-irp,Roip-iOj) and Ci^i2,iira,ri3, Rap) can be represented in the explicit 
form for the values of / given by conditions (3.46) and (3.47), respectively, (for details of 
calculation, see Appendix). For Tq, < Rap — r^, where < r^j < a^, we have 
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"^1 />-<^2 



X 3 ° f+,, } , / = /i + /2-2p>0, p=-l,0,l,...,p„„,, (3.52) 



r (/i + /2 + 1 



/ = /i + /2-2p+l>0, p = 0,l,...,p^a,, (3.53) 

where ji(x) = y^vr/ (2x)/;^i (x) and /i/(x) = y^Tr/ (2x)-K';_,_i (x) are the modified spherical 
Bessel functions of the first and third kind, respectively, [Q, T{z) is the gamma function, 

and ?/q,/3 = K-Ra/j. The dimensionless parameter t/q,/? = k-Rq/j [|?/q,/3| = 
v^-Ra/?/*^] characterizes the ratio of the distance between particles a and /? to the depth of 
penetration 5 of a plane transverse wave of frequency u into the fiuid. In a certain frequency 
range, the quantity Iva/sl can be both smaller (for short distances between particles) and 
greater (for space-apart particles) than one. 

Note that relations (3.52) and (3.53) for Fi-^i^ i^Va, r^j, Ra/3, i^) and Ci^i^^iiva, t^, -Rq/j) defin- 
ing the induced fiuid velocity and pressure were obtained without imposing any additional 
restrictions on the size of particles, distances between them, and the frequency range. 

In the particular case of equal spheres ( 1, 2, ... , N), if the observation 

point lies at the surface of a sphere (r^ = a) and rp = a, relation (3.52) for I = li + I2 — 2p, 
where p = 0, 1, . . . ,Pmax coincides with the corresponding relation (4.11) in pB| up to the 



factor (—1)^. For I = li + I2 + 2, relation (3.52) agrees with (4.11) in [|2B[ up to the factor 
{h + I2 + 5/2) in the first and second terms in (4.11) and the omitted factor 1/k\ (in terms of 



the notation used in [gg]) in the third term presented in (4.11) in [26| as the infinite series. 

For Va > Raf3 + ^/s, whcrc < rp < ap, we have 

2k {- ~ ~ TV'I'^ 

Tir]y 2 



X 



r(/2 + |)r(/i-/2-i 



a 



l = h + l2-2p>0, p= -1,0,1,..., p^ax, (3.54) 
22 



X 



r (/i - /2 - 1 



1 /^^aV 1 /^/^V 1 



2/l — 1 \Ral3 J + 3 \Ral3 J 2/ + 3 

/ = /i + /2-2p + l>0, p = Q,l,...,p^ax- (3.55) 



Note that the quantities F/iij.K^a, -Ra^s, t^) and Ci-^i^^ii^Va, rp, Rap) defined by (3.54) and 
(3.55) for Tq, > Rap-^-fp can also be represented in the form (3.52) and (3.53), which are vahd 
for Tq, < Rap — rp, performing the changes Tq, <-> Rap and h <-> I and putting, respectively, 
h = / + /2 - 2s, where s = -1, 0, 1, ... , s^^^; and s^ax = niin(^[(/ + /2)/2] , 1 + min(/, /2)), and 
Zi = Z + ^2 + 1 - 2s, where s = 0, 1, . . . , s^a^; and s^^j; = min(^[(/ + I2 + l)/2] , 1 +min(/, ^2)). 

The dimensionless parameter ap^ — ap/R^p (as well as Uap — Qa/Rap), which is the 
ratio of the radius of particles to the distance between them, is always smaller than one (for 
spheres of equal radii, it cannot be greater than 1/2). In diluted suspensions, apa <S 1. 

For I — > 0, the quantity Pi^2,i{''^a: Raf3: ^) defined by relation (3.42) with I2 — 2 
can be represented as follows (see Appendix): 

2 ~ 

Ph2,i{ra,ap,Rap,uj) = {-If — jh{xa)j2{bp)hi{yap) (3.56) 

TTT) 

for fa < Rap — cip and 

- 2 ~ 

Ph2,i{^a,ap,Rap,(^) = (-1)^+^ — hi^{xa)j2{bp)ji{yap) (3.57) 

TTT) 

for Va > Rap + dp- Here, bp = nap. 

In the particular case (3 — a, we have 

^K(^-r.,u;) = E T\Zl{ra,aa,u;) ■ fa^mM^ (3-58) 
pS(^-^«,a;) = E D\{Z{{ra.aa) ■ fa^mM^ (3.59) 

l2m2 



^Simi {Roc, = ^^^i2(ra, o^, o;) (^0„(a;) x e^^J , (3.60) 

<?ir(^- c.) = - /dr^e(a. - r^) ^ (r„, r^, c.) 
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• F^'"^>^\R^ + r:^,u;)Y*^M^^'a), (3-61) 

l2m2 

■ F^''>^>^\R^ + r'^,u;)Y*^M^<p'J, (3.62) 



where 



(3.63) 



2 f 

Fhh{ra,K,uj) = — dk 2 ^h{kra)3hikr^), (3.64) 

' 



(3.65) 



2 ^ 



Chi2ira,r^) = - dkkji^{kra)ji2ikr^), (3.66) 



TT 



oo 

2 /•„ k 



v^^^^'^^R^ + ro.^uj) = ^P,_2(ra,a„,cu)(0„(cu) x n<,), (3.68) 



Z f K 

Piih{ra,aa,i^) = — dkj^ — ^ jh{kra)jh{kaa)- (3.67) 

TTT] J K'^ -\- H'^ 
' 

Taking into account that according to (3.60) only one harmonic of i;^'^|*^^(i?a,ra,a;) is 
not equal to zero, we can represent this component as follows: 

9 

Q 

67r 

where Ua = 'Tajfa- 

According to (3.50) and (3.51), the quantities .f^;^mi,oo ^limi,oo nonzero only 

oo 

for at most three and two values of ^2, respectively. Therefore, the infinite sums Yl, in 

^2=0 

relations (3.58), (3.61) and (3.59), (3.62) are replaced by sums containing, respectively, at 
most three and two terms corresponding to different values of 12- This means that the 
quantities F;^;2(rQ,, r^, a;) and Ci^ijja.i'^'f^ should be determined only for these values of l-i- 
Finally, for Va > aa and < < a^, we can represent these quantities in the form (for 
details of calculation, see Appendix): 

Fi,Ur^,ri,u;) = (-1)^- \kMjU<) - (h - ^) 4 

l2^li + 2p>0, p = 0, ±1, (3.69) 

1 fr'V''^ 1 

Chhi'Tcx^r'^) = Si^,i^-i— [^] + Sr„,aSr^,a TTT i^hM+i " , /2 = ^1 ± 1 > 0, (3.70) 



'a \'ccj ^"a 
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where x ' = kt' 



a' 



For the quantity Pi^2{ra,cta,^), in the domain Tq, > a„, we obtain 

2 ~ 

PlA'^a, Cla, UJ) = hi{Xa)j2{ba)- (3.71) 

irrj 

In the particular (the observation point hes at the surface of sphere a) and 

= Qa, relations (3.69)-(3.71) are reduced to the form (see Appendix) 

F;,;,(a„,a„,a;) = (-l)''^j,_)(6„)/i;_(6„), h = h + 2p > 0, p = 0, ±1, (3.72) 

Cilia (««) «a) = ^ ((^ia.ii+l + '^«2A-l) ) ^2 = ^1 ± 1 > 0, (3.73) 

2 ~ 

Pl,2{aa, Cla, t^) = hi{ha)]2{ha), (3.74) 

Tirj 

where Z^aa: = max(/i,/2) and l^in = min(/i,/2). 

For I2 = h and I2 = ^1 + 2, expression (3.72) coincides with the corresponding result 
obtained in [^]. However, the quantity T^^^^H ) defined by relation (3.63) for = 

and = Oq is not equal to zero for three values of I2 [h = h, ^1 + 2, and /i — 2 (for li > 2)], 
while, according to p6|, Tj^™^( ) 7^ only for /2 = ^1 and I2 = h + 2. 

Relations (3.4), (3.6)-(3.8), (3.10), (3.12), (3.24), (3.26), (3.27), (3.29), (3.31), (3.32), 
(3.58), (3.60) and (3.61) completely determine the function v{r,uj) in the entire space pro- 
vided that the induced surface force densities are known. Analogously, relations (2.8), (3.5), 
(3.9), (3.11), (3.13), (3.15), (3.25), (3.28), (3.30), (3.33), (3.59), and (3.62) uniquely repro- 
duce the original function p{r, u) at any point of the space with the exception of the points 
of the surfaces of the spheres ( 1,2,..., N) where the obtained quantity is 

equal to the half-sum of the original function given at r = Ra + + and r = Ra + <2q, — 
because the original function defined as the fluid pressure for Tq, > and (2.17) for < 
is discontinuous at the surfaces = a^. With regard for (2.8), (2.17), and (3.5), the fluid 
pressure at the surfaces of the spheres can be represented as follows: 

p{Ra + flc. + 0, cu) = p^^°\Ra + aa, uj) + p'^^^iRa + a„ + 0, u;), (3.75) 

where 
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p'^^^iRa + a„ + 0, cj) = 2p'"'^(i2« + to) - iup {R^ + a^) ■ U^iuj) (3.76) 

and the quantities p^^\Ra + cia,uj) and p^"'^{Ra + aa,uj) are defined by the corresponding 
above-derived expressions for p^^\r,u) and p*"'^(r,u;) given at r = Ra + fla- 

IV. SYSTEM OF EQUATIONS FOR HARMONICS OF INDUCED SURFACE 

FORCE DENSITIES 

In the previous section, the required distributions of the velocity and pressure fields of 
the fluid as well as the forces and torques exerted by the fluid on the particles were expressed 
in terms of harmonics of the induced surface force densities. To determine these harmonics, 
we use the stick boundary conditions for the fluid velocity at the surfaces of the particles 

v{r,t) = Ua{t) + (n^{t) xr^ , a = l,2,...,N. (4.1) 

^ ^ r=Rc,+aa 

Passing in relations (4.1) to the Fourier transform with respect to the frequency and 
using representation (3.4) for the fluid velocity, we obtain 

V,(a„ u) = v'^^iR^ + a„, u;), « = 1, 2, . . . , iV, (4.2) 

where 

Va{aa,io) = Ua{uo) + (ri«(cj) X a^,) - v^^\Ro, + ao,,uj) (4.3) 



is the velocity of the point = Ra + o-a of the surface of sphere a relative to the fluid 
velocity v^'^^Ra + aa,uj) at this point in the absence of particles in the fluid. 

For each a, we expand the quantities in relations (4.2) in the spherical harmonics 
Yim{0,ip) analogously to expansion (3.19) for the induced surface force densities. Note that 
representations of the induced velocity and pressure of the fluid in the form (3.24)-(3.28) 
are, in fact, the expansions of these quantities in the spherical harmonics Yim{d,(f) at the 
surface of the sphere of an arbitrary radius Tq, with the center at the point Oa- Therefore, 
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putting Ta = tta in relations (3.24), (3.26), and (3.27), we immediately obtain the required 
expansion for the quantity v^"''^{Ra + o-a,^) at the surface of particle a in the spherical 
harmonics Yim{0,ip). As a result, we obtain the following system of algebraic equations in 
the unknown quantities f pi^{uj): 

N N 

E E T'JlZ (^) ■ fpMm. (^) = Vo^Mm. (^) - E {^Sl. (^) + (^) } ' (4-4) 

13=1 h-nii 13=1 

where 



{O)sol 



^ [dn^v^^>"\R^ + a^,uj)Y;„Xda,^a), (4.6) 



and, to simplify the representation of relations, we introduce the following notation: 



where the quantities Tll^^ia^^^ap, R^p^uj), Tj2™2(ao, a„, t^), 'y^*^"'^(i2„, a„, u;), 

(r)ind / -r-i \ (soViind / -r-, \ i (sol)ind / -w-t \ i j2 i 

<jimi(^«>«f.>^)> "^k/imi (^"'««'^)' and (i?a,a^,u;), are defined, respec- 

tively, by relations (3.34), (3.63), (3.31), (3.60), (3.32), and (3.61) for = Cq,. 

In the absence of the sum on the right-hand side of Eqs. (4.4) (this sum is absent in the 
case of the stationary problem and conservative external fields), the tensor T^'|^^^(u;) defines 
harmonic {hmi) of the expansion of the fluid velocity at the surface of particle a induced by 
harmonic {12^^12) of the expansion of the induced surface force density /^(a/j,^;) distributed 
over the surface of particle j3 (another particle if /? 7^ a or the same if /3 = a). For the 



stationary case (cu = 0), these quantities were introduced in ||3^ and called hydrodynamic 
interaction tensors. Their generalization to the nonstationary case was performed in p5| , |26| . 
In what follows, we use this terminology. Equations (4.4) agree with the corresponding 



equations (3.7) in ||2^ up to the terms v^^J^^^^ [uj) and v^[i°!IT'^ i^) ^ which is caused by 
extension (2.18) for the stress tensor mentioned above and the zero extension of the external 



forces to the domains occupied by the particles used in p6 . 
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V. STATIONARY CASE 



In the previous sections, we reduced the problem of determination of the velocity and 
pressure fields of the unbounded viscous fluid induced by an arbitrary number of spheres 
immersed in it to the solution of the infinite system of linear algebraic equations (4.4) in 
the harmonics of induced surface force densities fp^imi^)- this paper, we consider the 
important particular case corresponding to the stationary mode, i.e., all quantities are time 
independent. It is easy to perform the passage to this case using the general relations 
obtained above, setting a; = in them, and assuming that all quantities are independent 



of oj [for example, v'^"' {r,oo = 0) 



V 



ind/ 



r), T^^iu = 0) - T^^:^Z. etc., i.e., simply 



omitting the argument uj (or t)]. In view of relations (3.31) and (3.60), we get 



V 



(r)ind 
(3,1 



(R r ) = 

imi K-'-^ai I a) u. 



Mr 



0, 



1,2, 



{r)indf 



1,2, 



(5.1) 



(r) = 0, v'^^'^'^r) = v'P'''\r), (5 

Thus, in the stationary case, the fluid velocity caused by the induced volume forces is 
determined only by the solenoidal component of the external force field. This is quite natural 
because, according to (2.24), in the stationary mode. 



F^p'^\r) = -e(aa - ra)F(^°^)^^*(r). 



(5.2) 



Passing in relations (3.52), (3.54), (3.69), and (3.72) to the limit as a; —> 0, we obtain the 
corresponding relations for the quantities defining the fluid velocity in the stationary mode. 

For (5 ^ a, lim Fi^i^^i{ra,rfs,Rai3,uj), where < r/s < a/^ and I = h + k - 2p, p = 
-1,0, l,...,Pm„a,, gives that Fi^i^^i{ra,rfj, Rap) ^ only if Z = + ^2 (p = 0) and I = 
^1 + ^2 + 2 (p = -1) 



r (/i + /2 + 1 



^vRap r {h + i) r {k + i) ' 

Flil2,li+h+'2i^'a,rp, Rap) = Fhh,li+hi''"a,ri3, Rap) (ll +^2+2^ 



X{l-[li + l2 + - 



h + 



R. 



■aP, 



1/^/3 
^2 + 2 V R-ap , 



(5.3) 



(5.4) 
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in the domain < Rap — ff3, and only ii I = li — I2 {h > h, P = h) and I = h — I2 
2 (/i > /2 + 2,p = /2 + 1) 



TT 



/2 Tjh-h 



4r?r„ r (/2 + 1) r (/i - /2 + 



h > k, (5.5) 



TT 



4^vra r (/2 + 1) r (/i -h-f) 



X 



R. 



+ 



3 

/2 + ^ Vr„y 



h>h + 2 (5.6) 



_/l — /2 ~ I V Tq, / '■2 -|" 2 

in the domain Tq, > Raj^ + r^. 

Thus, in the stationary case, the sum over / in relation (3.34) contains only two terms 
with I = li + I2 and / = /i + Z2 + 2 for Tq, < R^p — and with I = li — I2 > Q and 
/ = /i — ^2 — 2 > for Tq, > Rap + rp. 

For /5 = a, for > fla, we have 

-F«i«2(^a5aa) = ^i2,h Fiiii{ra,aa) + Si2,ii-2Fij^i^{ra,aa), k = h + '2n>0, n = 0, ±1, (5.7) 
where 



Fl^h{ra,aa) 



Fi^,h-2{ra,aa) 



(2/1 + l)r]ra Vr„ 



1 



2r7ra Vr, 



/1-2 



h > 2. 



(5.8) 
(5.9) 



According to (5.7)-(5.9), 



FLio(aa,aa) = Si.i.FiJaa), l2 = li + 2n>0, n = 0,±l. 



(5.10) 



where 



-^Zi(aa) = Fij^ij^{aa,ao 



Therefore, the tensor 



(2/1 + l)r7a„ 



-iQ,Zim2 



(5.11) 



(5.12) 



is diagonal with respect to the indices li and I2 [0. Here, 
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-'■ a,hmi — (^2/1 + 1)^ -"-iimi.OO- {O.LO) 

We represent the system of equations (4.4) in the form 



o,/imi J a,;im2 a,limi ~r ^a,limi / y ) / y a,limi mi ( \" / 

m2=-h /S^^o I. '2^12 J 

separating terms corresponding to the hydro dynamic interaction between different particles 
on the right-hand side of the equations. In the absence of the external force field, system 
(5.14) coincides with the corresponding system derived in Taking into account the 

obtained explicit form for the quantities and Tf'j^™^, we get for any /i, ^2 > [ pOj] 



rpa,l2m2 r ^0 

This enables us to seek a solution of system (5.14) by the method of successive approxima- 
tions. As the zeroth iteration, we use system (5.14) with the zeroth right-hand side, which 
corresponds to the absence of interaction between the particles (the induced surface force 
density on the surface of particle a is determined by the characteristics of the fluid and this 
particle and independent of the characteristics of the rest particles). The solution of system 
(5.14) in the form of a series, each term of which, in fact, determines the contribution of a 
separate iteration to the total solution is given in Moreover, the final result is expressed 



in terms of the tensor T(^';^^^ and the inverse tensor T^^l^i defined by the condition 

h _ , 

^ a,h'mi ' ^ a,lim3 ~ '^mi,m2-'- 1^0. iOJ 

m3=—h 

Taking relation (5.13) into account, the inverse tensor T^^^^^ can be represented in the form 



^a,hm2 _ (2 /1 + l) ^a j^h 

3v/vr 



where -K'i^mi.oo tensor inverse to the tensor -K'jjmi.oo defined by the condition 

h ^ ^^^^ ^ 

-^hm-ifiO' ■^l\rn3,m~ ^rni_,m2l- (5.18) 



m3 = 
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However, the explicit form for the tensor T^l^^^ is not given in except for li = 0. 



Furthermore, it is stated that this tensor always exists due to the uniqueness of the solution 



of the Stokes equation. Being the central point for the determination of the solution in 



the derivation of the forces and torques exerted by the fluid on particles is based on this 
statement. As a result, these forces and torques are also expressed in terms of the inverse 

Ct li ?Ti2 

tensor T^\^^_^. For this reason, in the present paper, we dwell on the derivation of a similar 
solution in detail. To this end, we represent in the form 

/ a,limi / Q,Zimi) (5.19) 

n=0 

where f'^alimi a solution of system (5.14) corresponding the nth iteration. 

In a similar form, we represent the force and the torque Tq exerted by the fluid on 
particle a (as well as the induced velocity and pressure of the fluid) 

oo 

F^ = Y.F^:\ (5.20) 

n=0 

oo 

T„ = ^T("), (5.21) 

n=0 

where F^^^ and T*^"-* are the force and the torque corresponding the nth iteration 

F^:^ = -/So , (5.22) 
^1") = -^ E (e.x/SJ. (5.23) 



m=—l 



A. n = 0. Noninteracting Particles 

In this approximation (zero iteration), the inflnite system of equations (5.14) is reduced 
to the collection of independent systems of equations for each particle a and each li 

h 

m2=—h 

For h = 0, 
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-^00,00 ~ 3y^' — (5.25) 
and, according to (5.18) and (5.16), we have 

kZ,oo = 3v^, t:;;° = Cal. (5.26) 
We represent the solution faflo of system (5.24) with Zi = in the form 

f{0) _ Atfi) I j;{ea;t,0) /r ^yx 

/a,00 — /a,00 "T Ja,00 ) V^-^'J 

where the harmonics fafll and fa^o^^ associated, respectively, with the relative translational 
motion of particle a with the velocity 

[/W = Ua- v'""^ (5.28) 

and the external force field are determined as follows: 

/S = eaC/L^ (5.29) 
& - -ia^^^^t + ^r^"- (5.30) 

Here, v^^^^°'-{r)^°' — Vq^oo"' is the fluid velocity induced by the solenoidal component of the 
external force field in the absence of particles averaged over the surface of particle a and 

is the force acting by the solenoidal component i7'(*''')^^*(r) of the external force field F'^^*{r) 
on the fluid occupying the volume V^. 

Analogously to (5.27), in what follows, for any nth iteration, we represent the quantity 
fc^Jm a^ the superposition of the components 

j.(n) _ j.{t,n) n{r,n) p{ext,n) „ _ n 1 O (K 

Ja,lm — Ja,lm J a,lm J a,lm i — U, 1, Z, . . . yO.CiZ.) 

associated, respectively, with the translational motion of particle a, its rotation, and the 
external force field (as is shown in what follows, in the particular case / = 1, f'^a\m be 
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represented as a superposition of four components). To separate the contributions caused 
by the translational motion of particles, their rotation, and the external force field to the 
induced velocity and pressure of the fluid as well as to the forces and torques exerted by the 
fluid on particles, we also use analogous representations for these quantities. According to 
(5.27) and (5.32), for zero iteration, we have 

/S = . (5.33) 

Now we consider the case li = 1. According to (5.24), we have 

X! ^ImlflO ' /i.lma = { •\7I ^ ~ '"Mmi | ' (5.34) 

where 

j.sol _ Msol a{sol) /p. 

Using relation (3.48) for -K'ljmJ oo setting ^2 = = 1 in it, we obtain 

detK\Z,oo = 0- (5.36) 



Therefore, the tensors K\^^^ qq and T°|}^^ are degenerate and, hence, the inverse tensors 

1771(2 '~ ^ \tHj'2 

■^imi 00 ^a'lmi ^'^ ^'^^ cxist or they should be defined in another way |^ than by 
relations (5.18) and (5.16). Thus, even for the case of noninteracting particles, it is difficult 



to interpret the relation for obtained in |30| and expressed in terms of the nonexistent 



(in the ordinary sense) inverse tensor T^'^^^. Furthermore, with regard for the interaction 
between particles, not only harmonics with li = 1 but all harmonics ^^^^ are expressed in 



~a,lm2 



terms of T^'i^^ [relation (38) in ||3^]. To understand the reason leading to the impossibility 



of the existence of T^ -^^_^, we investigate the system of equations (5.34) in detail. To this 
end, first, we decompose all vectors in (5.34) into the independent unit vectors eo and e±i 
defined by relations (3.23) taking into account that any vector a = {ax,ay,az) may be 
represented in terms of these unit vectors in the form a = (a+i, a_i, oq), i.e., 

a = aoCo - (a^iCi + a+ie_i) , (5.37) 
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where a±i — {iay ± a^) / -\/2- 

Going from = (l]^^, fi^y, fi^^), 



(AO) N , sol ^ 

\J a,lm,,xj J a,lm,yj J a,lm,z J J cu-k-i "'a,lm — 



f(0) 



? / 

a,lm,a;' a,lm,j/' "'a,lm,^ 



U 



sol 



to _ (f^a,+l) ^a,-l) ^a,o)) — ifa,m+l->fa,m-l->fa,mo)-> ^^'^ 



c(0) 



f(0) 



c(0) 



? / 



sol 



If '^a mo)> reduce the system of nine equations (5.34) to three inde- 



pendent systems 



sol 



^Ja,00 Ja,+1+1 Ja,-1-1 — -'-'-'Sa"'a,f 



00 



/a,00 '^fa,+l+l + '^fa,-l-l — ^^a^afi + 10^q,M^°^,_i_|_i 



c(0) 



p(0) 



/a,00 + 2/(1^ 3/^ — iea^a,0 + 10Ca^a°-l-l' 



f(0) 



f(0) 



sol 



A AO) , AO) _ ■ ry 



/a,0-l + 4/(^J_io — ^£a^a,-l 10Ca'^a°+10 ' 



(0) 

«,-10 
p(0) 



' 0,0+1 



-10 ' 



where £q = 10aa^a/\/3, and 



(0) _ _12 f ,,30i 
Ja,+1-1 ~ 2 "Sa^aj+l-l) 

(0) _ _12 t ,,so« 



(5.38) 



(5.39) 



(5.40) 



(5.41) 



The solutions of systems (5.39) and (5.40) have the form 



u 



,sol 
a,-10 



AO) _ .^a^ 2 / 

/a,0+l — ~*"^"a,+l "T 2 ?a 

/a,0-l — ' g i'a,-l 3 l^"a,+10 ^"a,0-l J > 

■(0) - ,-£^0 , 4--^ (n,^"^ - Al,'"^ 
a,+10 — ^ g "q,-1 2 '5" l^a,0-l ^^a, 



,+10/ ' 



.(0) _ _ •£« Q I 2 t / soi _ 4 soi \ 



(5.42) 



The determinant of system (5.38) is equal to zero. Indeed, it is easy to see that the 
left-hand side of the first equation of system (5.38) is equal to the sum of the left-hand sides 
of the second and third equations of this system. For the consistency of system (5.38), the 
right-hand sides of the equations of this system must satisfy the following condition: 
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1 



Y: e„ ■ = . (5.43) 



m=— 1 



Using relations (3.61) and (4.6), we verify the validity of this condition because 

m=— 1 

E e„ ■ = • (5.45) 

m,=— 1 



Furthermore, relations (5.44) and (5.45) are also true in the general nonstationary case 
where u ^ including any values of and not only for = aa- 

Therefore, system (5.38) and, hence the original system (5.34), has an infinite number 
of solutions instead of a single solution as it is stated in . 



This result is quite natural. Indeed, the problem of determination of fa,im with the use 
of boundary conditions (4.1) means that we try to represent the unknown induced surface 
force densities /^(cia) in terms of the known fluid velocity at the surfaces of the particles. 
According to the continuity equation (2.2) for the incompressible fluid, the fluid velocity is a 
solenoidal vector. Thus, within the framework of this approach, the required quantity f ^{aa) 
can be determined only up to an arbitrary potential vector. The potential component of the 
induced surface force must make no contribution to the fluid velocity. Therefore, instead of 
the first equation in system (5.38), as an additional equation, we can use the condition of 
the absence of the potential component of the induced surface force density f^{aa)- 

In the general case where the interaction between the particles is taken into account, the 
conditions of the absence of the potential components of the induced surface forces can be 
written in the form 

Ya = 0, a = l,2,...,iV, (5.46) 



where 



— E ■ fa,lm " fa,00 + fa,+l+l + fa -1-1, a — 1,2, . . . , N. (5.47) 



m=— 1 
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In the approximation of noninter acting particles considered in this section, these condi- 
tions have the form 

= E • /SL = & + + fi%-i = 0, a = 1, 2, . . . , TV, (5.48) 

m=—l 

where the superscript in Y^'^\ just as in f^^jm other quantities, stands for zero ap- 
proximation. 

In the case where both the solenoidal and potential components of the induced surface 
forces should be determined, a certain additional equation linearly independent of the second 
and third equations of system (5.38) should be formulated. 

To derive additional equations with regard for the interaction between the particles in 
the fluid, we consider relation (3.5) at r = Ra and equate it to relation (2.17) given at 
r = Ra- As result, we get the relations 

Ya = AnVSal ^ {p^^^'iRa) + pf'^^'iRa)} , a = 1, 2, . . . , TV. (5.49) 

In the derivation of Eqs. (5.49), we used the exphcit form for the quantities £)|^™^(0, r^), 
where < aa, and obtained that p'^^^'^^{Ra) = 0. 

Putting Tq = £, where e — > -|-0, in relations (3.25), (3.28), (3.38), and (3.40), we can 
represent the quantity p~p^^'^{Ra) in the form (for details of calculation of the required 
quantities for = 0, see Appendix) 

-If oo h ^2+1 

I2=lm2=-l2 m=-{l2+l) 
XYi2+i,miQap,^ap)^, (5-50) 

where n«/3 = Rap/Ra/s- 

Thus, in the general case, Eqs. (5.49) differ from Eqs. (5.46) corresponding to the ab- 
sence of the potential components of the induced surface forces by the nonzero right-hand 
side. However, in the approximation of noninteracting particles, equating the right-hand of 
Eqs. (5.49) to zero, we obtain Eqs. (5.48), which means that the induced surface forces have 
no potential components in this approximation. 
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Using Eq. (5.49) and the second and third equations of system (5.38), we obtain the 
system of three equations with nonzero determinant, the solution of which has the form 



[•(0) _ ,-^"r> t ( A^.SOI I „.soi 



sol 



fa -1-1 — —i-^^afi ~ 3 [^'^'^a-l-l + ^a,+l+lj ) 
/a,00 — o ^" V"a,H 



+1+1 "T -1-1 I • 



(5.51) 



Returning in relations (5.41), (5.42), and (5.51) to the quantities fa\mi get 



AO) ^ f{r,0) n 
J a,lm J a,lm ' J 



a.lm 5 



(5.52) 



where 



J a,lm — 2 a ^ ' 



/{extfi) 
a.,\m 



-^ea(4<L+(<L)^), 



(5.53) 
(5.54) 



and = S^J'^c^Q is the Stokes friction coefficient for a rotating sphere of radius a^. The 
quantities [u^a^im) defined as follows: We consider three vectors bm = {hm,x-, bm,y, bm,z): 
where m = 0, ±1. According to (5.37), we go to the vectors bm = {bm+i,bm-i,brao) and 
represent bm in the form 



where 



1 

m2=— 1 



^00 —bo-i —bo+i ^ 



mi = 0, ±1, 



(5.55) 



B 



m\m2 



^+10 ~^+l-l ~^+l+l 
6_10 -&-1-1 -&-1+1 

Then bZ, is the vector defined as follows 



V 



mi, m2 = 0, ±1. 



(5.56) 



mi=— 1 



mi = 0, ±1, 



where B^^^^ is the matrix transposed to the matrix 
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(5.57) 



By virtue of (5.52), the translational motion of particles makes no contribution to the 



harmonic [ f^'^ 

In view of relation (5.54) for fa^m \ a similar representation may be given for any right- 
hand side of system (5.34) for which this system is consistent. Indeed, we can analogously 
represent the quantity fa',im 



0" 



J a,lm 



(5.58) 



where the vector b^"^^ = (fia x e^) can be rewritten in the form (5.55) with the asymmetric 



tensor B^;^^^^^ 



a,mim2 















—i^lafi 










mi, 1712 = 0, ±1. 



(5.59) 



According to (5.59), we have (sJLim,)^ = -^£Um., which gives (b'^^^Y = -b'^^l,. 
For li > 2, according to (5.24), we have 



(5.60) 



where 



j.sol ^ Msol , a{sol) 
a.lm aim ~ a.lm 



(5.61) 



~ hm2 



Assuming the existence of the inverse matrix -f^^j^^^oo defined by relation (5.18) for 
/i > 2, we obtain 



^(0) _ Aextfl) ^ _ j,a,hm2 ^^g^l 

m2=—l\ 



mi a,/im2' 



(5.62) 



i.e., in the approximation of noninteracting particles, the harmonics of the induced surface 
force densities faj^mi > 2 are not equal to zero only if the external force field has the 

solenoidal component. 

Substituting the obtained solutions (5.27), (5.29), (5.30), (5.52)-(5.54), and (5.62) into 
(3.58), (3.59), (3.29), and (3.30) and using (3.34), (3.35), and (3.43)-(3.51), we obtain the 
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following relations for the t, r, and ext components of the harmonics of the fluid velocity 
and pressure caused by translational motion of particles, their rotation, and the external 
(solenoidal) force field, respectively, in the approximation of noninteracting particles: 



(5,t,0)m<i/„ N 
*^a,hmi {-"-a, 'a) 



V 



{S,r,0)ind. 



{S,extfi)ind / tj 



{S,t,Q)ind, r> 
Pa,hmi {-"-o 



{S,ext,0)ind , 



Pa,li 



mi 



Rr 



V 



{S,t,0)ind 



{Ra 



{S,rfi)indi „ 



(S,ext,0)ind / tj 



{S,t,0)ind/ Tf 



{S,ext,0)ind , 



Rn 



(S,r,0)indf p 
P/3,hmi 



1 

iv 10 



Kr 



Jjit) 



E'T'hm2( \ Aext,0) 

7712= — h 

1 ^ X \ X "^Q * i-r(t) 



1712 = 



Ta) = CpKmiira, ap, R^p) ■ uf, /3 ^ ft. 



2 a/3 



E a/3, i^a/j) • (o^ X e;;, J , 

/3 7^ a. 



m2=— 1 



/3,/2m2' 



l2m2 



l2m2 



(5.63) 
(5.64) 
(5.65) 
(5.66) 

(5.67) 

(5.68) 
(5.69) 

(5.70) 

(5.71) 
(5.72) 

(5.73) 



,ra)=0, /3 = 1,2,...,7V. 

According to (5.73), the rotation of the particles has no effect on the fluid pressure 

'^'°^''^V) - 0, /3 = 1,2,...,7V. (5.74) 

Substituting (5.63), (5.64), and (5.66) into (3.24) and (3.25) and taking relation (3.76) 
into account, we obtain the following space distributions for the corresponding components 
of the velocity and pressure fields of the fluid: 



da ^ 




1 - 


(^)l 











(J-3n„n„)^C/W, 



vi^'^-'^'^'^'^iR^ + r«) 



p(f'*'°)-^(il. + r.) = ^%n..t/W. 



(5.75) 

(5.76) 
(5.77) 
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Relations (5.75) and (5.77) coincide with the well-known distributions of the fluid velocity 
and pressure induced by a sphere moving with the constant velocity [/„ in the immovable 
fluid (t;*"-^ = 0) or with the velocity and pressure flelds of the fluid uniformly moving with 
the velocity v*"-^ relative to an immovable sphere {Ua = 0) ||^. Relation (5.76) describes 
the known distribution of the velocity of the fluid induced by a sphere rotating in it with 
the constant angular velocity 0^ p. 

Using relations (5.68), (5.69), and (5.71), we obtain the following main contributions to 
the fluid velocity and pressure induced by particle j3 in the vicinity of the surface of particle 



a for Va <^ Rap (the near zone): 

^(sm^nd^j^^ + r„) ^ ^a/3a (I + Uapriap) ■ Uf, (5.78) 

^is,rfl)tnd^j^^ + r«) apal^ {Uf^ x riap) , (5.79) 

j>f'^'^-\Ra + ra)^lv^n^,-Uf. (5.80) 



However, for the correct description of the fluid velocity and pressure up to these powers 
of the parameter a, these quantities should be also found in the flrst approximation. 

Using relations (5.5), (5.6), and (3.55), we obtain the velocity and pressure flelds of the 
fluid induced by particle /? 7^ a in the far zone, i.e., for Tq, ^ Rap- In zero approximation 
with respect to the ratio Rap/ra, we obtain the following space distributions of the fluid 
velocity and pressure induced by a system of spheres far from it: 

o N 

^(5,t,oW(^^ + r J ^ (/ + n.nj ■ ^ a^uf, (5.81) 

^'^a 13=1 

v^'^'^-'^'^'^^Ra + ra)^-Y.4 X n„) , (5.82) 

^a 13=1 
q 1 N 

pism^nd^Ji^ + r.) ^ E «/3 ^« ■ Uf. (5.83) 

^ 13=1 

In this approximation, can be interpreted as the distance from a certain point inside this 
system (e.g., its center) to the point of observation. 

Naturally, these relations can be also derived by using the superposition principle for the 
velocity and pressure flelds of the fluid induced by noninteracting particles keeping only the 
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main terms in them. 

In the case of equal particles ap — a, relations (5.81)-(5.83) can be represented as follows: 
^ismind^R^ + r„) P. |- (/ + n«n«) • C/*°*, (5.84) 



a 

3 



^(5,r,0)md^_^^ + re,) ^ (fX*''* X U^) , (5.85) 



' a 



p(sm^nd^ji^ + r„) « ^4 n« • t7*°*, (5.86) 

2 ri 



where 



N 

J7*°* = ^[/W, (5.87) 

13=1 

N 

= ^ np. (5.88) 

/3=1 

Thus, the velocity and pressure fields of the fluid induced by the system of equal spheres 
in the far zone, in zero approximation with respect to the ratio of the typical distance 
between two spheres to the distance to the point of observation, can represented as the 
velocity and pressure flelds of the fluid induced by a single sphere moving in the fluid with 
the translational velocity [/*°* + v^"'^ and rotating with the angular velocity f2*"*. In the 
particular case of two equal spheres rotating in the opposite directions with equal angular 
velocities, their rotation has no effect on the fluid velocity far from the spheres within the 
framework of the considered approximation. 

If all spheres move with the same translational velocity Up — Uq {U^p = t/^*-* = 
C/o -v'"-'^), then 

^^smind^R^ + r,.) « (/ + n^n^) ■ C/^, (5.89) 

4 To, 

j,is,mrnd^R^ + r.) 5^ • [/W, (5.90) 

Z r 

N 

where atot — J2 ci/s- In this case, the action of the system of spheres in the far zone is 

/3=1 

equivalent to the action of a single sphere of radius atot- 

If all spheres rotate with the same angular velocity = Q,, then we obtain the following 
relation for the fluid velocity: 
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^(5,r,ow^j2^ + r„) ^ (n X n„) , (5.91) 

Air 

which corresponds to the fluid velocity induced by a single sphere occupying the volume 

N N 

Vtot = E = (47r/3) E 4 far from it. 

/3=1 /3=1 

According to (5.22), (5.27), (5.29), (5.20), and (5.33), in the approximation of noninter- 
acting particles, we obtain the following relations for the components of the force acting by 
the fluid on particle a caused by translational and rotational motions of particles and the 
external force field: 

= -aC/i*\ (5.92) 
= 0, (5.93) 

Relation (5.92) is the classical Stokes force acting by the fluid on a sphere moving in 
it with the constant velocity Ua if the fluid moves with constant velocity d*"-^. Relation 
(5.93) illustrates the well-known fact of the absence of a force acting on a sphere due to 
its rotation. The first term in relation (5.94) corresponds to the known Faxen relation |1| 
defined by the classical Stokes law with the velocity equal to the velocity of inhomogeneous 
motion of the fluid in the absence of spheres averaged over the surface of sphere a. 

Substituting (5.92)-(5.94) into (2.27) and taking into account the absence of the inertial 
force -F^ defined by (2.33) in the stationary case, we obtain the following relation for the 
total force acting on particle a: 

ptot ^ pe.t _ ^ ^(,^0) ^ ^^^iO)sol^^f"^ (5.95) 

where 

~ext ~(p)ext ~{sol)ext f , . . , , 

=jdrF^^\r) (5.96) 

Vac 

is the force acting by the external force field F'^^^{r) on the fluid occupying the volume V^. 
The difference F'^^^ — F^ in (5.95) is the force acting on particle a in an arbitrary force 
field F^^^ir) corrected for the buoyancy force for this force field. 
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To determine the torque T^' exerted by the fluid on sphere a, we substitute relations 
(5.52)-(5.54) into (5.23) and perform certain transformations taking into account the validity 
of the following relation: 

1 1 
J2 {em xbl) = - J2 i^rn X b„) (5.97) 

m=—l m=—l 

for any vectors hm and 6^ defined by relations (5.55) and (5.57). As a result, we get 

Ti*'°) = 0, (5.98) 
T(:'°) = -ef fl., (5.99) 
Ti-*'°) = 2a(a„x^(0M(r))''" - ti^°'^'^^*. (5.100) 

According to (5.98), we have the well-known result of the absence of the torque acting on a 
sphere moving with a constant velocity in an unbounded fluid. Relation (5.99) is the classical 
result for the torque exerted by the fluid on a sphere rotating in it with constant angular 
velocity Q,a [|l[. The first term in (5.100) corresponds to the torque exerted by the fluid on 
sphere a due to the inhomogeneous motion of the fluid in the absence of particles caused by 
the solenoidal component of the external force field. For a homogeneous distribution of the 
fluid velocity, this torque is equal to zero. The second term in (5.100) defined as follows: 

is the torque exerted on the fluid sphere occupying the volume due to the solenoidal 
component F^''°'-^''^\r) of the external force field F'^'^^r). 

Substituting (5.98)-(5.100) into (2.28), we obtain the following relation for the total 
torque exerted on sphere a: 

Tf = - f"^' + T(:'°) + 2i^(a^ X ^;(o)-'(r))^^ (5.102) 

where 

tr = f^t" + T^r'^'" = Jdr (r„ X F-(r)) (5.103) 
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is the torque acting due to the external force field F^^*(r) on the fiuid occupying the volume 
Va- The difference T^^* — T^^* in (5.102) is the torque acting on sphere a in an arbitrary 
force field F^^^{r) corrected for the buoyancy torque for this force field. 

Thus, the external force field makes the contribution to the total force and torque exerted 
on sphere a both due to the buoyancy force and torque on the one hand, and due 
to the inhomogeneous distribution of the fiuid velocity in this force field in the absence of 
particles, on the other. 



B. n = 1. The First Iteration 

For the first iteration, with regard for (5.24), the system of equations (5.14) is reduced 
to the form 

h ( h ] 

Erpa,hm2 _ /;(!) — _ J rpl3,l2m2 _ ^(0) 0{sol) I 

m2=—h I3^cx ym2=—l2 ) 

where the quantities fa]2m2 determined by relations (5.27), (5.29), (5.30) for I2 — 0, 
(5.52)-(5.54) for k = 1, and (5.62) for k > 2. 

Taking into account relations (3.34), (3.35), (3.37), (5.3), and (5.4), we represent the 
solution of system (5.104) for Zi = in the form (5.32) with n — 1, where 

= E ^pTMiRap) ■ Uf, (5.105) 

/So = E Cp(^p(^a/3(Tpa (^a/? X ftp) , (5.106) 

S,3' 



p(ea:t,0) 1 
I3,hfn2 J ' 



oo h 



+ T. E (5.107) 



l2=l m2 =—l2 



where 



TM{Rap) = Tl'?^{aa, a/3, Rap) — o 5 1 + '^aP'^afi) 

+ + Crla) - riapriap^ | (5.108) 
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is the modified Oseen tensor . 

According to (5.105)-(5.107), in this approximation, faflo ~ o"^, fa'ol contains terms 
proportional to a and a^, and oo^'' terms ~ 1/ EJ^n, where = 1, 2, . . ., furthermore, 



For li = 1, system (5.104) has the form 



(5.109) 



I^lm2 -f(l) — \^ J rpl3,l2m2 AO) I „,/3(soO I I'Cilin'l 

/ . -"-Imi.OO J a,lm2 ~ rz 2^ \ a,lmi ' J I3,l2m2 ^ " a,lmi ( i \O.LL\J) 



m2=-l V ^ fj^a \^m2=-l2 

which differs from analogous system (5.34) for the iteration of noninteracting particles only 
by the right-hand side. This implies that system (5.110) is consistent only if 

E ' 51 ] E ^a,lmi^ ■ fp}2m.2 + '^a!lml (5.111) 

mi=-l I3j^a \^m2=-l2 ) 

This condition is satisfied because the following relations are true: 

E e„,-TS^=0, (5.112) 

mi=— 1 

E e^,-<S^=0. (5.113) 



,/3(so«) 

mi=— 1 



Therefore, system (5.110) has an infinite number of solutions that differ only in the 
potential component. To determine the unique solution, it is necessary to introduce an 
additional equation linearly independent of the equations of this system. To derive this 
equation, we use relations (5.47)-(5.50) and get 

E e^./S^ = 4vrv^a^E{pr"'(^")+#°^^"'(^^^ a = 1, 2, . . . , iV, 

m=-l /3j^a 

(5.114) 

where pf'^^'"'^{Ra) is defined by relation (5.50) with /^^/^^^ ffj2m2- 

Unlike Eq. (5.48), Eq. (114) contains the nonzero right-hand side. Therefore, the inter- 
action between the particles leads to the appearance of the potential components f^^\a) of 
the induced surface forces densities /^(a). 
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We solve the system of equation (5.110), (5.114) in a similar way as in the case of 
noninter acting particles. The final result has the form 

fW _ f{tA) I f{r,l) I Aext,l) ,0,1) f'S 115"! 

J a,lm J a,lm J a,lm ' J a,lm ' J a,lrm \o.i.LO) 



J a,lm 






Ap,i) 

J a,lm 


y(i) 

— f,* a 
Cm 3 , 




"a, mi 


\ " \ " rpl3,l2m2 _ 
/ y / y o,lmi 

fSj^a l2m2 


,(t,0) _ 
J I3,l2m2 




\ " \ " rp/3,l2m2 _ 

/ ^ / ^ a,lmi 


fir,0) 

J f3,l2Tn2 



where 

I , e — t,r, ext, (5.116) 

(5.117) 

T.^pTi'Z^-Uf, (5.118) 

rp/3,l2m2 _ j;(r,0) 
a,lmi ■ J /3,l2m2 
jSj^a l2m2 

= I: E E X W?^,,2m) >^2m(e.;3, ^a/s), (5.119) 
"^'l ISyta m=-2 

T(ext,l) _ j rp0,l2i-n2 Aextfl) B{sol) I /r , r)r.\ 

"a,mi ~ S a,lmi ' J p, 121112 ^a,lmi r ) {O.IZU) 

Pyia (l2m2 } 

and the quantities (ba',m) determined by relations (5.56), (5.57), and (5.118)-(5.120). 

In the first iteration, we obtain that /^^'i^ ~ cr^, has terms proportional to cr^ and 
(7^, and fa^il^^ contains terms ~ l/R^^, where n — 2,3, 

For h > 2, the solution of the system of equations (5.104) can be represented in the form 
(5.32) with n = 1, where 

_ _ V- e V- rjp»'hm2 0^00 jAt) Cc;i91^ 



I3^a m2 = -h 



V3 V- ' 



,(r,l) ^ _ V£ ^Q,/im3 _ „/3,lm2 /q x P* 

J a,limi n / J / J / J o.,l\m\ a,hm3 \ P m2 

^ ^+a. m3=-/i m2=-l 

= -v^(2/,+i)Ef ^i^v,„ I: E k;;:;,oo 

/37^a m2=-h m=-(;i+l) 

• X r,,+i,„(e„^, $„,,«), (5.122) 

/3^am3=-/i I ^ i O m2=-l 



/3^a m3=— ii 

a^>> oo ^2 I 

.,(ext,l) /,(ext,l)y I _ rr:P,l2m2 f(extfl) _ l3{sol) I /C^IOQ^ 

^"/3,m2 "T- J f a,hm3 ' J /3,l2m2 " a,hms [ ■ \0.i.^0) 

(2=2 m2=— (2 
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Thus, the interaction of particles leads to the appearance of all harmonics for the com- 
ponents of the induced surface force density connected with the translational motion of 
particles and their rotation. 

Since the potential component f^'^\aa) of the induced surface force density f^^\aa) 
has only harmonics with I — 1 defined by (5.117), we obtain the following results for the 
corresponding fluid velocity and pressure {va > Oa) induced by these potential forces: 

^{S,p,i)ind^^^ = 0, /3 = 1, 2, . . . , TV, (5.124) 

pf'^'^'^^'ir) = ^'^^^"'^(iia + a„ + 0), /? = 1, 2, . . . , TV, (5.125) 

where 

+ a, + 0) = --^5^. (5.126) 
47rV3a^ 

Thus, the potential component f^'^\ap), where P — 1,2,. . .N, of the induced surface 
force density makes no contribution to the fluid velocity, while the fluid pressure caused 
by this potential force is equal to zero everywhere with the exception of the surface of the 
sphere where this force is distributed. 

In view of (5.117), the potential component f^^'^\aa) is represented in the form 

n VW 

f^''^ = T^r^, (5.127) 

and, hence, has only the radial constant component. In this case, it is easy to see that the 
potential component of the induced surface force has no effect on the forces and torques 
exerted by the fluid on the particles immersed in it. 

These conclusions concerning the influence of the potential component of the induced 
surface force densities on the fluid velocity and pressure as well as on the forces and torques 
exerted by the fluid on particles also remain valid for any nth iteration if the potential 
component /^'"^(a^) of the induced surface force density f^^\aa) corresponding to this 
iteration has only harmonics with I — 1 defined by a relation similar to (5.117). 

By virtue of relations (5.105), (5.106), (5.115), (5.116), (5.118), (5.119), (5.121), and 
(5.122), in this approximation, the main contribution to the fluid velocity and pressure 
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in the vicinity of particle a due to the interaction between the particles caused by the 
motion of particles and their rotation is given by the quantities D(^^'*'i)*"'^(r), t;(^5,r,i)indj^^^^ 
p<^,t,i)tnd^j,^^ g^j^j^ p(5,r,i)jndj-y.j^ Hioreover, 7;(^^'*'i)*"'^(r) and t)(^^''^''i)*"'^(r") are of the same order 
in the parameter a that the quantities 'y^^'*'°^*"''^(r) and v^p'^''^^^"'^{r), where /3 7^ a, defined, 
respectively, by relations (5.78) and (5.79) in the approximation of noninter acting particles. 
Near the surface of particle a, for <C Rai3, the main terms for the fluid velocity and 
pressure due to the hydrodynamic interaction between the particles can be written as follows: 

+ n^pn^p)-Uf, (5.128) 
^ vf^'''^'-\R^ + r.) + vf^'^^'-^R^ + r.) «^ - -) (/ - 3n„n«) • ^ H 

X (7|, (n«^ X ftp) , (5.129) 
E pf'^^'^^'iRa + r-a) +p(f'*'^)^"''(i2. + r.) p(f'*'^)^"'^(H« + r J - (l - 1) 

^ E ^/3a {I + riapn^fs) ■ Uf, (5.130) 

^(5,.,iW(^^ + r«) = (1 - J) 1^ n« • E «/3^§a {n.p x fl^,) . (5.131) 

According to (5.128)-(5.131), the main terms of the fluid velocity and pressure caused by 
the hydrodynamic interaction between the particles due to their translational motion and 
rotation are proportional to the flrst and second powers of the parameter cr, respectively. 

Substituting (5.105) and (5.106) into (5.22) and (5.116), (5.118), and (5.119) into (5.23) 
and taking relation (5.97) into account, we obtain the following relations for the forces and 
torques exerted by the fluid on sphere a for this iteration corresponding to the translational 
motion of particles and their rotation: 

F^'^^ = E ipTuiR^p) ■ Uf, (5.132) 

F^'^^ = - E C/3a/3(7a/3(7/3a (^a/? X ^p) , (5.133) 



8r, 
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T'^a'^ = -iaaa E ^"/3^/3- >< Uf) , (5.134) 

j.ir,i) = _k ^ (/ _ Sn^f^n^p) . Qf,. (5.135) 



Relation (5.132) agrees with the second term in relation (44) in corresponding to 



the first iteration. It is worth noting that relation (5.133) for the force exerted by the fluid 
on a sphere due to rotation of the rest spheres coincides with the known result for two 



spheres |T^, while the corresponding force given by the first term in relation (49) in ||3y] is 
proportional to the inverse tensor -fiT^^^^ qq, which, as was mentioned above, does not exist. 
Relations (5.134) and (5.135) also agree with the known relations for the torques exerted by 
the fluid on a sphere due to the translational motion and rotation of the rest spheres |]l],|l2|. 

According to (5.132) and (5.108), the force F^^'^^ contains terms proportional to the first 
and third powers of the dimensionless parameter a. Restricting ourselves only to the terms 
proportional to a, we obtain 

F'i''^ ~ E (I + riapn^p) ■ Uf. (5.136) 

The terms proportional to in (5.132) should be retained only if the problem is solved 
within the framework of at least the third iteration. 

Often, it is necessary to determine the velocities of particles in a fluid in given external 
force fields (in particular, the sedimentation velocity of particles in a fluid due to the action 
of the gravity force) and to investigate the influence of hydrodynamic interactions between 
the particles on their motion. Let [/„ and fla be, respectively, the velocity of translational 
motion and the angular velocity of a single sphere a in a fluid caused by the action of a 
certain external force field. In the case of several spheres, these velocities change due to 
hydrodynamic interactions between the spheres. Using relations (5.133)-(5.136), in this 
approximation, we can represent the changed velocities Ua and fl^ as follows: 

Ua = Ua + Y.uil (5.137) 
= 0„ + E (5.138) 
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where U^J^p and Q,^^^ are, respectively, the changes in the velocity of translational motion 
and the angular velocity of sphere a due to the motion and rotation of sphere /?. Up to the 
main terms in powers of the parameter a, these quantities have the form 

t^S = cr/3a 1^ (i" + no^prio^p) ■ Uf - afsattp {riap X ^^^)| , (5.139) 
= 'l^Pa 1^ K/3 X uf) + al^ - n^pn^p^ ■ fl^j . (5.140) 

The first term in relation (5.139) agrees with the main term in the relation given in |jl| 
for the problem of sedimentation of two particles with constant velocities without rotation 
under the action of the gravity force. The first term in relation (5.140) coincides with the 
well-known result in [|l| for the angular velocity of sphere a, whose rotation is induced by 
the motion of sphere (3 moving in the fluid with the relative translational velocity uf. 

C. n = 2. The Second Iteration 

Using Eqs. (5.14), (5.49), and (5.50), we obtain the following systems of equations for 
the nth iteration, where n > 2: 

h 

Erpa,h'm2 _ An) _ _ /j-i/3,«2m2 . An-l) /r i ^i \ 

-'■ a,hmi ' J a,lim2 2-^ 2-^ a,l\m\ ' J /3,/2m2' yo.LHiLj 

m2=—h I3y^al2m2 

Yl""^ = E . = ^^VSal E P?'"-'^^"'(i?.), (5.142) 

m=—l P^a 

where pf''^~^^^"''^{Ra) is defined by relation (5.50) with ff^{^ substituted for ff^i^^. 

Putting n = 2 in (5.141) and (5.142), we obtain the required equations for the second 
iteration. Further, putting /i = and using relations (5.105)-(5.107) obtained for the first 
iteration, we obtain the following solution: 

Py^a [ 7^/3 m2=-l 

■ Hz + -EE T^;,r ■ /&I2 1 > (5-143) 

^ /2=2m2 = -/2 j 
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/a,'00 = X! 1 TMiRali) ' ^ ^7^7^/37^7/3 (^7 ^ ^/37) + o C/3 X! ^a^OcT^ 
/3y^a I 7y^/3 ^ 1712=-! 



( rp^ CO I'Z I 

• + } - E E T^^r ■ f';:lL > (5-144) 

Aext,2) _ f rjiP,hm2 /?(ext,l) /j. , 

/a,00 — ?a 2_/ a,00 ' J (3,l2m2- {O.L^O) 

P^a l2m2 

According to (5.143) and (5.144), the quantities fafll and f^^^ are infinite power series 
in the parameter a starting from the second and third powers, respectively, which essentially 
differ them from the corresponding relations (5.105) and (5.106) obtained for the first itera- 

(t 2) 

tion. The quantity f a',oo contains only even powers of a. We note that the representation for 
a harmonic of the induced surface force density as an infinite power series in the parameter 
(7 is also vahd both for other I (Z = 3, 4, . . .) and for higher iterations (n = 3, 4, . . .). 
For li = 1, the system of equations (5.141) has the form 

1 d h 

ET^lm2 J!(2) _ rpl3,l2m2 n{l) /p. ^.(.\ 

^ ImifiO ' J a,lm2 — f^^ 2^ 2^ a,lmi ' J j3,l2m2-> 
m2=-l V /3^am2=-i2 

which is similar to system (5.110). Since the determinant of this system is equal to zero and, 
by virtue of relation (5.112), this system is consistent, it has an infinite number of solutions 
that differ only in the potential component. To obtain the unique solution, we add to this 
system Eq. (5.142) with n — 2. Analogously to the first iteration, we represent the solution 
in the form 

x(2) ^ At,2) Ar,2) Aext,2) »(p,2) /p. 



where 



/&! = e;^, (5.149) 
^>Sa = E E T'Xl ■ ^ = i^r, ext. (5.150) 

P^OC l2Tn2 

According to (5.148) and (5.150), the quantities fa'Jmi and fa'Jln are infinite power 
series in the parameter a starting from the third and forth powers, respectively. 

51 



For li > 2, the solution of system (5.141) can be represented in the form (5.32) with 
n = 2, where 

Ae,2) _ _ V- V- f^'hms f},l^m2 . .(6,1) . _ ^ ^ ^5 1 5ll 

J a,h-mi — «,iimi a,hm3 J/3,«2m2' t — i,f,fiXi. {O.LOl) 

P^a 1713=— h l2m2 

The quantities faf^mi ^^"^ fa'hmi infinite power series in the parameter a starting 
from powers of (2 + Zi) and (3 + /i), respectively. 

Analogously, we can determine the harmonics of the induced surface force density (and, 
hence, the velocity and pressure fields of the fluid) for higher iterations. 

In the investigation of analogous problems, the main question is connected with the 
substantiation of the reduction of the infinite system of equations in unknown quantities 



27| , |28[| (for the considered approach, the unknown harmonics faim)- Within the framework 



of this approach, we can reformulate this question as follows: How many iterations must 
be carried out to obtain the velocity and pressure fields of the fluid, the forces and torques 
exerted by the fluid on particles immersed in it up to a^, where p is a certain positive integer? 
Using the results obtained above, we can show that for any iteration n > 1, the main terms 
of the quantities fa',hn^ f^aTrn \ f^'afm are proportional, respectively, to cr'+" and o-'+"+-^, 
moreover, starting from n = 2 [for fa^^\ starting from n = 0], these quantities contain 
infinite sums of higher powers of a. For this reason, to obtain the induced velocity and 
pressure fields of the fluid up to a^, it is necessary: 

(i) to carry out p iterations, 

(ii) for each sth iteration, where s < p, to retain only harmonics with I = 0,1, . . . ,p — s, 

(iii) in each retained harmonic, to retain all terms up to terms proportional to inclusively. 
Further simplification of the results depends on the point of observation. 

Substituting (5.143) and (5.144) into (5.22) and (5.147), (5.148), and (5.150) into (5.23) 
and taking relation (5.97) into account, after certain transformations, we obtain the following 
relations for the forces and torques exerted by the fluid on particles due to their translational 
motion and rotation corresponding to the second iteration: 
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Pjta I 77^/3 m2=-l 

• + (CJl -EE ^^Sor • f';£j, (5.152) 

I 7//? m2=-l 



• + } + E E i^Sr-zSiU k (5-153) 

Tt''^ = E C/J^X^ 4 E ^7 [n^P X (TM(i«/37) • + I ^/^- E E f^Klm 

/3^Q 17^/3 V ) ^ m2=-lm=-2\ 

X (46(*i), + (6(*iy ) y.„(e.„ $.,) - ^ E E E < 

/ i2=2TO2=-/2 m=-(Z2 + l) 

X X ll2+l,n.(e.;3, ^a/3)}, (5-154) 

^a'^^ = -«a E '^a/?! ^0.^0.^^ otfi(y fio. {J " n^fi^a.^) ■ ^ot " ©(-/V " 3) E ^7«7^/37^7/3 

P,»n.2 



O 1 2 / rr 



X (n^p X (n,;, X fi,)) - -e^ja^a E E J W^o^2m x {^h%^, + (bj'^'^) 



m2=— 1 m=— 2 
00 ^2 '2+1 

xy2^(e«;3,$a/3) + 47rE4 E E (i^&i,™ >< /KU 

^2=2 m2=-l2 m=-(i2+l) 

X Yi2+i,^(0„^,<l'„;3)|. (5.155) 

If we restrict our consideration to the second iteration, then only the main terms in 
powers of a should be taken into account in relations (5.152)-(5.154). As a result, we get 

+ e{N - 3) (I + Ua-yUa^) ■ E ^/37 (^ + np^^fs^) ■ uf >, (5.156) 

F^a ^ ~ "^07 j ^aO'a^a'^^ia {'^af X fi^) 

+ e{N - 3) {I + na-yUa^) ■ E ^/3«/3^/37^7/3 (^/37 X f^/3) k (5.157) 

/3^7,o J 

77^a V 
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{iI + np^n^^)-U^^)\\, (5.158) 



7^0 



/3^7,« J 



(5.159) 

The first terms in (5.156)-(5.159) correspond to the self-interaction of particle a due to 
the action of the fluid induced by this particle and reflected from the rest particles. The 
second terms in these relations, which are nonequal to zero only for > 3, correspond to 
three-particle interaction. For fixed a and /3, these terms describe the contribution to the 
force and torque acting on particle a by the fluid induced by particle (3 and scattered by 
all rest particles (except for a). For a system of two spheres, relations (5.156)-(5.159) agree 
with the well-known results given in p|JT^ 

Taking relations (5.132)-(5.135) and (5.156)-(5.159) into account, we can represent the 
velocity of translational motion Ua and angular velocity JIq, of particle a in a fluid in given 
force fields in terms of the velocities Up and the angular velocities fi^ of noninteracting 
particles immersed in the fluid in these fields as follows: 

N 

Ua = Ua + Y.Uap, (5.160) 

13=1 
N 

0=1 



For (3 ^ a 



3 f 3 







7^a,/3 

+ apap^a^p {I + n^^rio,^) ■ {np^ x ^0)^, (5.162) 

fla0 = ^^i^j + 6(A^ - 3) ^a7^/37^7a|a/30-/3T,(n„^ X [Up^ X fi^)) 

+ ^ {(^^7 X + (^/37 ■ (^"7 X ■^/37)}| (5.163) 

are, respectively, the velocity of the translational motion of particle a and its angular veloc- 
ity induced due to motion and rotation of particle (3. The first terms in (5.162) and (5.163) 
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characterizing two-particle interaction are determined by relations (5.139) and (5.140), re- 
spectively, describing the corresponding changes in the velocities of particles within the 
framework of the first iteration. The second terms in relations (5.162) and (5.163) charac- 
terize the changes in the translational velocity of particles and their angular velocities due to 
three-particle interaction. The power of the parameter a contained in these terms is greater 
than that in the first terms by one. The second iteration not only changes the quantities 
Uai3 and f2Q,/3 but also leads to the appearance of new terms Uaa and Q,aa equal to 

3 f 3 1 

Uaa = -T E ^^7^7" I 4 '^-'^ + 3n„^n„^) ■ C/^*) + a^a^^ (n„^ x Q,^) \ , (5.164) 

3 r 3 1 

^ca = -. E '^a7^7a ot^f ^-^{^1 X (n^^ X fi^)) + - (n^^ X U^^^ \ , (5.165) 

which are absent for the first iteration and characterize the changes in the velocity of transla- 
tional motion Ua and the angular velocity Vta of particle a due to its motion in the presence 
of other (not necessary moving or rotating) particles. The second term in (5.165) agrees 
with the expression given in [|T| for the angular velocity of a sphere that can freely rotate 
in the fluid due to its motion with constant velocity Ua in the presence of other spheres. 
Note that relations (5.164) and (5.165) also follows from the terms in relations (5.162) and 
(5.163) proportional to 0(A^ — 3) by formally replacing G(A^ — 3) by 1 and putting (3 = a. 



D. Friction and Mobility Tensors 

Within the framework of the second iteration, using relations (5.22), (5.23), (5.92), (5.93), 
(5.98), (5.99), (5.132)-(5.135), and (5.156)-(5.159), we represent the forces and F^^^ 
and the torques T^*-* and T^^'-* exerted by the fluid on particle a as follows: 

(5.166) 
(5.167) 
(5.168) 



a 


= E F^^''^^ = 

k=0 


N 

-Ea- 

/3=1 


^ (3 


a 


= E F^-'^ = 

k=0 


N 

~ E ^ap 

(3=1 




rp{t) 

a 


= E n^'^ = 

k=0 


N 
/3=1 


Tj{t) 

^ (3 
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2 N 

T'^^ = E ^i^''^ = - E C • (5-169) 

fe=0 13=1 



3 

Cla = -CaOa J E ^a7'^7« " n^j) , (5-172) 



where and are, respectively, the translational and rotational friction tensors and 
^Q,^ and are the friction tensors that couple translational motion of particles and their 
rotation. These quantities are infinite power series in the dimensionless parameter a. Up to 
the main terms corresponding to the second iteration, these quantities have the form 

^IJ = ^aM + AlJ), (5.170) 

C' = eaM + Af«), (5.171) 
3 

— ~^aa — (^la ) > (5.173) 

+ e{N J2 ^ocjCrjpaf}^{{e-nf)^)-naj{na-yXn^^)}\, P ^ a, (5.174) 

-e{N-3)j (^ajcrjacrp'y{{e ■ riaj) + {np^ X Ua^) np^}\ , P ^ a, (5.175) 

where e is the absolutely antisymmetric unit tensor of the third rank, {e ■ n) — eijkrik, 

9 

^11 ^ 16 ^ f^aT^^Ta + ^najU^j) , (5.176) 

= I E <^7a (J - ria^ria^) , (5.177) 



-Q{N-3)j (^7a^/37(-f + ^a7»^a7) ■ (-^ + ^/37'^/37)[' (5.178) 

+ e{N-3)^ Y (—) afs^a^pal^Uria^ ■ np^) I - ni3^na^\\. (5.179) 
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Taking the explicit form of relations (5.170)-(5.179) into account, we can show that the 
friction tensors satisfy the Onsager symmetry relations |T| 



{e,iy=eai (c^)^=c {^?iy=ej, a,(3=i,2,...,N, (s.iso) 

1 means the transposition of the matrix with respect to the space variables 
hj = x,y,z. 

Relations for the translational friction tensors ^^J, where P = 1,2, . . . , N, coincide with 
the corresponding well-known relations given in pOf retaining in them terms up to the second 
order in a inclusively. 

According to (5.170)-(5.179), the friction tensors of translational (^^J) and rotational 
i^ap ) motions of particles are determined up to and cr"^, respectively, and the tensors 
and are determined up to a . In the particular case of two particles, the friction 
tensors defined by (5.170)-(5.175) with N = 2 calculated up to the above-mentioned orders 



in a agree with the results given in [O] 



In view of (5.166)-(5.169), the force and the torque exerted by the fluid on particle a 
due to the translational motion and rotation of all particles are defined as follows: 

N 



= + = - E {ej . Uf + ej ■ O,} , (5.181) 

/3=1 
N 

r„ = T« + tW = - E • + ej ■ ■ (5.182) 



0=1 



Solving the system of equations (5.181), (5.182) for the quantities uf and O^, we obtain 



N 



= = - E ■ + f^l^ ■ r,} , (5.183) 

13=1 
N 

13=1 

where the translational (/^^J) and rotational (fJ-^^) mobility tensors as well as the tensors 
/x^^ and n^J that couple translational and rotational motions of particles determined up 
to the same orders in the parameter a that the corresponding friction tensors have the form 
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If 3 1 

A*IJ = 7- + (1 - Safi) (I + ■n^pnaf3)j , (5.185) 



1 r n-3 



/^«" = ^ l'^-/^^ + (1 - ^ (3n./3n,;3 - /)| , (5.186) 



= (1 - 5a/3) (e ■ n,^) , (5.187) 

/^fj = mI^'- (5.188) 

It is easy to see that the mobihty tensors defined by relations (5.185)-(5.188) satisfy the 
Onsager symmetry relations |jl| 

(/xJJ)^ = mIJ, (/xjff = /xf?, {f^^^f = til^, a,/3 = l,2,...,iV. (5.189) 

Within the framework of the considered approximation, for which ^t^J, Ai^/f , MaJ? ^^id 
fi^^ are determined up to terms of cr^, a'*, a^, and o"^, respectively, relations (5.185)-(5.188) 
for the mobility tensors coincide with known results given in . Unlike the friction tensors 



calculated in this approximation, the mobility tensors are determined only by two-particle 
interactions and fi^^ = At^J = 0. In order to take into account three-particle interactions 
in the mobility tensors, it is necessary to carry out calculations for higher iterations. 



VI. CONCLUSIONS 

In the present paper, we have proposed the procedure for the determination of the time- 
dependent velocity and pressure fields of an unbounded incompressible viscous fiuid in an 
external force field induced by an arbitrary number of spheres moving and rotating in it 
as well as the forces and torques exerted by the fiuid on the particles. The corresponding 
quantities are expressed in terms of the induced surface force densities. We showed that the 
relations for the required harmonics of the induced surface force densities given in I^O] for 
the stationary case are expressed in terms on nonexistent inverse tensors. We analyzed in 
detail the reasons for zero of the determinant of the corresponding system of equations. We 
formulated the consistent system of algebraic equations in the harmonics of these induced 
surface force densities. In the stationary case, we obtained relations for the harmonics up to 
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the second approximation inclusively. The obtained general results for the fluid velocity and 
pressure and the friction and mobility tensors corresponding to the second approximation 
agree with the well-known results obtained by other methods in various particular cases. 

The proposed procedure can be used for the investigation of hydrodynamic interactions 
of particles in the fluid for higher-order approximations in a similar way as it is realized in 
the present paper. At the same time, the relations for the fluid velocity and pressure and the 
forces and torques exerted by the fluid on the particles expressed in the present paper in terms 
of induced surface forces and the proposed procedure for the determination of these forces 
can be regarded as a basis for the study of hydrodynamic interactions between particles 
in the nonstationary case. Certain results concerning the time-dependent hydrodynamic 
interactions of particles in a nonstationary fluid will be given in subsequent papers. 
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To determine the explicit form of the quantities Fi^i^^i{ra-i rp, Rai3,uj), Ci^i2,i{ra, r^, Rap), 
and Pi^i.^^i{ra,ai3, Ra/3,uj) defined by relations (3.36), (3.40), and (3.42), it is necessary to 
take integrals of the product of three spherical Bessel functions. For this purpose, first, we 
consider the integral 



where a is real, i^, and 7 are real quantities nonequal to negative integers. Re k > 0, and 
a, b, and c are real positive quantities such that c > a + b. The case where the order of a 
Bessel function is a negative integer is reduced to the above-considered integral using the 



APPENDIX: 



00 



,a-l 




(Al) 
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known relation J^n{x) = (— valid for integer n. To take this integral, we use the 
Hankel method |35|. To this end, we consider the following auxiliary integral: 

+ 



fdx -^^J^{ax)J^{bx) \h^'\cx) + {-iyH^'\-cx)] , (A2) 





where H!^^\x) is the Bessel function of the third kind and p is a certain quantity. For 
p = l — a — u — fi, integral (A2) can be reduced to the form 



OO 1 



dx— -J,{ax)J^{hx)H^^\cx). (A3) 



We consider the integral 



dz—--Uaz)J^{hz)H^^\cz). (A4) 
c 

in the plane of the complex variable z along the closed contour C consisting of the large Cr 
and small Cr semicircles of radii R and r, respectively, centered at the origin of coordinates 
and lying above the real axis, and the segments [— -R, — r] and [r, E\ along the real axis as 
i? — > 00 and r — > 0. Taking into account that the functions Jv{z) and HI^^z) are analytic 
functions of z in the entire complex plane of z with the cut (—00, 0], we obtain that integral 
(A4) is determined by the simple pole of the integrand at the point z = in provided that 
the integrals along the contours Cr and Cr are finite. Since R 00, for c > a + 6, by the 
Jordan lemma, the integral along the contour Cr is equal to zero if a < 4^. In the case 
where c = a + b, the integral along the contour Cr is equal to zero for a < 3| and finite for 
a = 3|. The condition of finiteness of the integral along the contour as r ^ imposes 
the following constraint on the quantity a for 7 7^ 0: 

a > 7 — z/ — yU. (A5) 

Furthermore, the integral along the infinitely small semicircle Cr is nonzero only in the case 
of equality (A5). In the particular case 7 = 0, inequality (A5) is strict. 

Passing in (A4) to the limit as i? 00 and r — 0, we can represent integral (A2) in the 
form 

60 



X 



a-1 



dx-^ -JJax)Ja{hx) H^^Hcx) + i-lYH^^X-cx) = 2(iK)'''W+^"^ UaK)lJbK)KJcK) 

+ L ' 'J 



p = l — a — z/ — /i, a > '-y — fi — u, 



(1 + z), (A6) 



a < if c > a + 6 or a < Oc if c = a + 6, 



where 



4i, if c> a + 6, 



Or 



3|, if c = a + 6 



(A7) 



and lyix) and Ky{x) are the modified Bessel functions of the first and second kinds, respec- 
tively. 

To determine the required integral (Al), we note that in (A2) 

"/f(^Hcx) + {-\YIl'^^\-CX)\ = [l + (_l)7---M-«] J^(cx) + [l + (_l)7---A*-«+l] Y^icx), 

(A8) 

where 1^1, (x) is the Bessel function of the second kind. In the particular case where a = 
7 — z/ — /i + 2k, where k = 0, 1, 2, . . ., condition (A5) is true and we obtain the following 
result for the required integral (Al): 



dx 



X 



J^{ax) Jf^{bx) J.y{cx) = (m)" ^^j^+m ^j^(^aK)I^{bK)K^{cK) + Sa^-y-u-fi (1 ~ (^7,0) 



2"-^r(7) a^b" (-1)^ 



(A9) 



a = 7 — z/ — /x + 2A;, /c = 0, 1, 2, . . . , 

a < if c > a + 6 or a < ac if c = a + 6. 

In the particular case where j = i', c> a + b, and a = 2k — fi < 9/2, where k = 1,2, . . ., 
relation (A9) agrees with the known result given in p. 232]. 

Setting in (A9) u = li + ^, n = I2 + \, and 7 = / + |, where li, I2, and / are nonnegative 
integers, we get 
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f X ~ ~ — 3 



a-3 



X 



Til + \ 



+ ^c,a+b ^aA 



r {h + 1) r (/2 + 1) 

a = / - /i - /2 + 2A;, A; = 0,1, 2 



8a6c 



, ^, . . . , 



(AlO) 



a < as if c > a + 6 or a < if c = a + 6, 



where ji{x) = y7r/(2x)/;_^i (x) and = y it/{2x)Ki_^_i{x) are the modified spherical 

Bessel functions of the first and third kind, respectively, |Q and 



«3 = < 



5, ii c > a + b, 
4, if c = a + 6. 

Passing in (AlO) to the limit — > 0, we get 



(All) 



(ixx""-^ ji,{ax)ji^{bx)ji{cx) = 7r22'~('i+'2+3) 



r(/i + |)r(/2 + | 



2 V2/i + 3 



+ 



\~ Oc,a+b Oa,4: 



■l)'=7r 



(A12) 



2^2 + 3 2/2 -ly J c'+i ■ Sabc 
c> a + b, a,b> 0, a = I - h - I2 + 2k, = 0, 1, 2, . . . , 

a < 0:3 if c > a + 6 or a < if c = a + 6, 

except for the special case where a, h, I2, 1 = 0. 

For c > a + b, relation (A12) agrees with the result given in pO, p. 239]. 

Note that we can apply the method used above for the determination of integral (Al) to 
the integral 



dxx" ji^{ax)ji^{bx)ji{cx). 



(A13) 



In a similar way as for integral (Al), we get 



dx ji^ {ax)ji^ {bx)ji (ex) = vr 2 2 



a— 5 



r(/ + i) 



r(/i + |)r(/2 + |) c'+i ■ """-^"""^ 8abc 
c> a + b, a,b > 0, a = I — li — I2 + 2k, k = 1,2, . . . 



a < as if a + b or a < if c = a + 6. (A14) 
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Relation (A14) coincides with (A12) for > 1. At the same time, we see that the 
direct apphcation of this method to (A13) does not enable us to determine integral (A13) 
for a = I — li — I2 (fc = 0) because relation (A14) is true for A; > 0. 

Analogously, we can determine the integral of the product of two spherical Bessel func- 
tions. The final result is presented as follows: 



00 







x'^ + K 





'1 + 








'2 + 


i) 





6>c>0, a = li-l2 + 2k, fc = 0,1,2,..., 
a < a2 if h > c or a < a2 if b = c, (A15) 



where 



4, if 6 > c. 



a2 = { (A16) 
3, if 6 = c. 



For b > c and > 0, integral (A15) agrees with the result given in p6| , p. 213]. 



Passing in (A15) to the limit as k — 0, we obtain the known Weber-Schafheitlin integral 
°? I2 Y (l-i + -] ( 1 / /)2 \ 1 



b -^[^2^2 



+ V'5«,3^^^, (A17) 
6>c>0, a = /i-/2 + 2A;, A; = 0, 1, 2, . . . , 
a < a2 if 6 > c or a < a2 if b = c, 



except for the special case a, li, I2 = 0. 



Relation (A17) agrees with the results given in |3^, pp. 239 and 209] obtained for b > c 
and b = c, a = 3, respectively. 

Note that except for the particular case where b = c and a < 3, it is simpler to determine 
the integrals of the product of two spherical Bessel functions (A15) and (A17) on the basis of 
relations (AlO) and (A12) for the integrals of the product of three spherical Bessel functions 
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found for a,b,c> 0. For this purpose, it is necessary to pass to the hmit a = (or 6 = 0) in 
(AlO) and (A12) defined for a < and to take into account that ji{0) = Si^o and the fact 
that for a = (or 6 = 0), the condition of the finiteness of integral ( A4) along the contour 
Cr leads to the decrease in the maximum possible value of a by one, i.e., as a2- After 
the corresponding renaming, we obtain, respectively, (A15) and (A17) defined for a < a2- 

In the particular case 6 = c, it is convenient to represent integrals (A15) and (A17) as 
follows: 



X 7 + <Ja,3 775 — ' 



a-1 



r (j'max ~^ 
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— Ifnin Imax + 2/c, A) — 0, 1, 2, ... , (Al8) 
+ 5^,3 , a ^ Imin - Imax + 2k < 3, k^ 0,1,2,..., (Al9) 

where Imax = max(/i, k) and Imin = min(/i, I2). 

To determine the quantity Fi^^ij^^ra, r^, Rai3, defined by relation (3.36) for I — I1 + I2 — 
2p > 0, where p — —1, 0, 1, . . . ,Pmax, we use the above obtained integrals (AlO) (for rp ^ 0) 
and (A15) (for = 0) setting a = 2 in them. Indeed, in this case, the representation of 
a in the form a = / — /i — /2 + 2k, where A; = 0, 1, 2, . . ., corresponds to A; = 1 + p. Since 
> > 0/3, depending on the distance to the point of observation, two cases are possible, 
namely, r„ < — and Tq, > R^fi + rp [the domain i?^/? + '"/j > > — should be 
considered separately because it does not satisfy the condition c > a + 6 of integral (AlO)]. 
Depending on the case, R^^^ or should be taken as the quantity c in relation (AlO). As a 
result, we obtain 

2k{~ ~ ~ 7r^/2 r (/i + /2 + f 



2 r (/i + Ij r (I2 + 



X 3",;^,, , l^h + l2-2p>0, p^ -1,0,1,..., Pmax: (A20) 
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for ra < Ra/3 - rp and 

2k. (~ ~ ~ TT^/^ 

irr) y z 



r(/2 + |)rf/i-/2-i) xy^-^ 



27 V 2, 
l = h + l2-2p>0, -1, 0,1,..., Pmax, (A21) 

for ra > Rap + rp. Here, Xa = KVa, xp KVp, and 7/^/3 = /«-Ra/3- 

Note that in view of the representation (AlO) for the integral of the product of three 
spherical Bessel functions, relation (A21) for the quantity Fi^i^^i{ra,rp, Rap,oj) valid for 
fa > Rap + fp can also be represented in the form (A20) valid for Ta < Rap — fp with 
the changes Tq <-> Rap and li <-> I putting li — I + I2 — 2s (i.e., p — I2 — s), where 
s = -1, 0, 1, . . . , s^ax and s^ax = min([(/ + l2)/2] , 1 + min(/, k) 

Putting rp = and Tq, = Rap in (A20) and (A21) and equating the obtained relations to 
one another, we obtain the following relation for the modified spherical Bessel functions: 

ji{x)hi+2{x) - ji+2ix)hi{x) = ^ + ^) ' ^ ^ 0- (A22) 

Note that this relation can also be derived by equating relations (A. 15) for c = 6 and 
(A. 18) putting a = 2 and li = I2 + 2 in them. 

To determine the quantity Ci^i^^i{ra-iT p. Rap) defined by relation (3.40) for I = li + I2 — 
2p+ 1 > 0, where p = 0, 1, . . . ,Pmax, we use relation (A12) with a — 3 and take into account 
that rp < ap. As a result, for rp ^ 0, we get 

3 1 rhr^2 



^ r /i + /2 + f) r[lr. 



Q.Ura,rp,Rap) - ^,,+,+1 {k + l)r {I2 + 1) R'lt^' ^^''^ 



for ra < Rap - rp and 



Q.,2Ara,rp,Rap) - -Yr{i2 + i) ^"^^^ I r(;,-;. + |) + ^'''^"'^-^ 



X 



r (/i - /2 - i 



1 \ 1 rp \ 



2k-l\Rap 2l2 + 3\Rap 21 + 3 
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(A24) 



for ra > Ra/3 + rp. 

Just as for the representation of the quantity Fi^i^^iira^rj^^Rap^oj) for > Rap + ?^/3, 
we can show the possibihty of the representation of relation (A24) for the quantity 
Ci^i2,i{ra-iTp,Rap) in the form (A23) vahd for < Rap — with the changes Tq, ^ Rap 
and li ^ I putting, li — i + ^2 + 1 — 2s, where s — 0, 1, ... , s^ax ^-nd s^ax — 
min([(/ + /2 + l)/2] , 1 + min(/, h)). 

In the particular case — 0, to determine the quantity Cij^i2,i{'ra, 0, Rap), it is necessary 
to use relation (A17). As a result, we obtain 

Qii2,/(ra, 0, Rap) = Si,^oCi,,i{ra, Rap), / = ± 1 > 0, (A25) 

where 

Cl,,l{ra, Rap) = i^Sl,h+l + ^r^,Rc.0 \ i^hh-l " ^l,h+l)^ (^26) 

for Ta < Rap s-nd 



G„Kr„, Rap) = ^ Um-1 {^f ' + Sr^,R„, \ iSi,i,+, - Sl,l,.^) \ (A27) 



a 



for ra > Rap- In this case, the quantities Cij^,i{ra, Rap) should be considered only for 
/ = /i ± 1 > because for I2 = 0, Pmax = if /i = or Pmax = 1 if /i > 1. 

Note that simply passing to the limit in (A23) and (A24) as — > 0, we obtain 

Chi,,i{ra, 0, Rap) = Si,,oCi,,i{ra, Rap), l^h±l>0, (A28) 

where 

CUra, Rap) = Si,i,+, ^ (^] (A29) 

for ra < Rap and 

Ci„i{ra, Rap) = Si,i,.,^ ( ^ ) (A30) 

for ra > Rap- 
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^aP \Rocp, 



1 / Rap ^ '1 



Comparing relations (A29) and (A30) with relations (A26) and (A27), we see that the 
function Cii^i{ra, Ra/s) corresponds to the first term of the function Ci^^iira, Rap)- Further- 
more, in view of (A29) and (A30), the function C/^,/(ra, Rap) has a discontinuity at the point 
Ta — Rap-i while, according to (A26) and (A27), the function Ci^,i{ra, Rap) is continuous at 
this point and equal to 

Ci,,i{Ra(3, Rap) = + , Z = /l ± 1 > 0. (A31) 

In a similar way, using (AlO), we can find the quantity Pii2,i{T^a,0'P, Rap,^) defined by 
relation (3.42) with I2 — 2. Indeed, according to (3.31) and (3.51), this quantity should be 
determined only for / = /i ± 1 > 0, which agrees with the constraint for the quantities /i, I2, 
I, a, and k in the integral defined by relation (AlO). As a result, we get 

2 ~ 

Ph2,i{ra,ap,Rap,i^) = (-1)^ — jh{^a)j2{bp)hi{yap) (A32) 

TTT] 

for Ta < Rap — dp and 

2 ~ 

Ph2,i{r-a,a0,Rap,uj) = (-1)^+^ — hi^{xa)j2{bp)ji{yap) (A33) 

TTTj 

for > Rap + ap. Here, hp — nap. 

Using the obtained integrals of the product of two spherical Bessel functions (A15) 
and (A17), we can determine the quantities Fi^i^{ra-iT'^,uj), Cij^h{i^a,i^a), and Pi,2(?'a, Oa, <^), 
defined, respectively, by relations (3.64), (3.66), and (3.67), for > and < < a^, 
taking into account that /^(^a? ''"q) and Ci-^i^{ra,r^) should be determined only for 
I2 — h + > 0, where p = 0, ±1, and ^2 = ± 1 > 0, respectively. As a result, we get 

T^V ( V 2/ x% \ra^ 

/2 = /i + 2p>0, p = 0,±l, (A34) 
1 fr^V''^ 1 

Cl^l^ira, r^) = (5;2,;i_i— I — + Sr^^a Sr^,a TT^ {^l2,li+l " <^/2,/i-l) ) 

Z2 = /l ± 1 > 0, (A35) 

2 ~ 

Pl,2{^a,aa,'^) = hi{Xa)j2{ba). (A36) 

TTT) 
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In the important particular case where = and = a^, using relation (A22), we 
reduce relation (A34) to the form 

Fi,i,{a^, a«,u;) = (-1)^ — j;_(6a)/i;_(&a), ^2 = /i + 2p > 0, p = 0, ±1. (A37) 

We also need the explicit form for the quantities Ci-^i^^i{ra, ap, Rap) and Ci^i^^i{ra, cia) for 
= 0. Setting Tq, = in (3.40) and using (A17), we obtain 

Ci,i,,i{Q, ap, Rap) = Si,,o / = /2 ± 1 > 0. (A38) 

Note that this result can be also obtained from the general relation (A23) for 
Ci^i2,i{ra-i Tpi Rap) found for > if we set rp — ap in it and pass to the limit as — > 0. 

To determine the quantity Ci^i^{ra-i cia) defined by relation (3.66), where I2 — h^l > 0, 
for Ta = 0, we, first, find the quantity Ci-^^i^{ra-i ««) for < tta- To this end, using relation 
(A17), we get 

Ci,i,{ra, a„) = 6i,,i,+^ \{—)'\ ^2 = ± 1 > 0. (A39) 



Passing in (A39) to the hmit as — > 0, we obtain 



1 



Y\-m Ci^i^{roc, Oa) = (^ii,o ^h,! -r- {MQ) 



At the same time, according to relation (3.66) for the quantity ^1/2(^0, ^'a), we obtain 
the following relation for = 0, = a^, and ^2 = ± 1 > 0: 

^1^2(0,0^) = 5i^fl5i^,i sinxl^^^) (A41) 

that has not any definite limit. 
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